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Introduction 


For  improved  care  of  patients  with  acute  lung  injury,  existing  clinical  data  sets  are 
examined  to  characterize  quantitatively  the  pulmonary  pressure-volume  (p-V)  curves.  A 
mechanistic  model  of  the  total  respiratory  system  (TRS)  is  constructed  for  both  the  in¬ 
flation  and  an  ideal  deflation  process.  The  mechanistic  model  computationally  simulates 
changes  in  the  TRS  based  on  an  application  of  statistical  mechanics  to  a  very  large  number 
of  elements  comprising  the  TRS.  Based  on  the  model,  various  shapes  of  the  p-V  curves 
are  related  to  the  corresponding  changes  in  intrarespiratory  conditions,  as  well  as  to  the 
magnitudes  of  p-V  curve  parameters. 
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Body 


Research  contents  are  reported  in  four  chapters: 

Chapter  1.  Quasi-Static  Pressure- Volume  Curve:  Comprehensive  Data  Analysis 
reports  results  corresponding  to 

Objective  1  (Examination  of  accuracy  and  limitations  of  the  sigmoidal  equation), 
Objective  2  (Development  of  a  method  for  quantitative  characterization  of  p-V  curves), 
in  Statement  of  Work  of  Research  Proposal. 

Chapter  2.  Mechanistic  Model:  Part  I.  Model  Development  for  Inflation  Process 
reports 

derivation  of  the  mechanistic  model  for  the  inflation  corresponding  to 
Objective  3  (Development  <5f  a  mechanistic  respiratory  model) 
in  Statement  of  Work  of  Research  Proposal. 

Chapter  3.  Mechanistic  Model:  Part  II.  Examination  of  Clinical  Data 
reports 

examinations  of  p-V  inflation  curves  based  on  the  mechanistic  model  as  part  of 
Objective  2  and  3 

in  Statement  of  Work  of  Research  Proposal. 

Chapter  4.  Mechanistic  Model:  Part  III.  Deflation  Process 
reports 

development  of  mechanistic  model  of  an  ideal  deflation  process  comparisons  with  the 
measured  data  sets  as  part  of 
Objective  2  and  3 

in  Statement  of  Work  of  Research  Proposal. 
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Chapter  1.  Quasi-Static  Pressure- Volume  Curve: 

Comprehensive  Data  Analysis 


Abstract 

A  p-V  model  equation  with  four  parameters  is  used  to  represent  various  existing  (p-V) 
curves.  The  report  is  focused  on  the  case  in  which  the  equation  is  applied  to  two  existing 
groups  of  p-V  data  (one,  twenty  nine  p-V  curves  of  healthy  adults  and  the  other,  twenty 
one  p-V  curves  of  patients  with  acute  respiratory  distress  syndrome)  to  determine  the 
magnitudes  of  the  parameters  for  each  data  set.  The  equation  is  found  to  represent  the  p- 
V  curves  of  both  data  groups  extremely  well.  It  is  also  confirmed  that  the  magnitudes  of  the 
four  parameters  of  the  error  function  p-V  model  equation,  combined  with  the  corresponding 
normalized  representation  of  p-V  curves,  quantitatively  distinguish  different  respiratory 
conditions  between  the  two  groups  as  well  as  between  different  data  sets  in  each  group. 
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Nomenclature 


pressure  (interpleural  pressure  difference) 

(volume-)  gradient  pressure  range,  Eq.(3) 
lower  (upper)  corner  pressure,  Eq.(4) 
pressure  at  maximum  compliance  increase  (decrease),  Eq. (5) 
pressure  at  the  inflection  point  (at  the  maximum  local  compliance) 
where  V  =  {Vv  +  VL)/2 
non-dimensional  pressure,  p/po  —  1 
volume 

lower  (upper)  volume  asymptote  (Fig.l) 

non-dimensional  volume,  (V  -  ((Vu  +  VL)/2)/(AVj2),  (Eq.(lb)) 
Vv  ~  Vl  (Fig.l) 

Greek  symbols: 

constant  defined  in  Eq.(la) 
apoAV (non-dimensional  parameter)  (Eq.(lb)) 

Ap/2  (Eq.(lb)) 


P 

Vgrad. 

Pcl(u) 

Pmci(d) 

Po 

P 

V 

VHU) 

V 

AV 


a 

A 


u 


Acronyms: 


ARDS 

LIP 

TRS 

UIP 


acute  respiratory  distress  syndrome 
lower  inflection  point,  Eq.(6) 
total  respiratory  system 
upper  inflection  point,  Eq.(6) 
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Introduction 


Quasi-static  pulmonary  p—V  (pressure  -  volume)  curves  provide  quantitative  informa¬ 
tion  on  the  respiratory  system  that  is  important  for  both  research  and  clinical  guidances. 
A  typical  inflation  p-V  curve,  obtained  for  an  anesthetized  human  subject  in  supine  posi¬ 
tion,  consists  of  a  nearly  linear  region  of  high  compliance  (i.e.  large  dV/dp )  sandwiched 
between  two  segments  with  low  compliance  at  low  and  high  pressure  regions.  The  shape 
of  the  curve  is  affected  by  two  mechanisms,  the  distension  of  the  elastic  respiratory  wall 
tissue  components  and  the  recruitment  of  the  alveoli  (‘pop-open’  mechanism).  The  latter 
is  the  opening  of  alveoli  overcoming  the  surface  tension  at  the  interface  between  the  gas 
and  the  liquid  film  lining  the  alveolar  surface.  A  pressure  increase  (i.e.  an  increase  in  the 
interpleural  pressure  difference)- results  in  the  recruitment  of  a  greater  number  of  alveoli. 
The  high  compliance  is  believed  to  be  associated  with  both  the  distension  of  open  parts 
and  the  (alveolar)  recruitment  of  collapsed  parts  of  the  total  respiratory  system  (TRS). 

In  order  to  quantify  the  characteristics  of  p  —  V  curves  as  well  as  their  changes  ob¬ 
served  in  clinical  settings,  various  p-V  model  equations  have  been  proposed  [1-8].  One 
commonly  used  model  equation  is  developed  by  dividing  the  entire  p-V  curve  into  three 
regions,  a  high-pressure,  low-compliance  upper  region,  a  high-compliance  midregion  and 
a  low-pressure,  low-compliance  lower  region.  The  midregion  is  represented  by  a  linear 
equation  between  p  and  V;  while,  the  two  low-compliance  regions  are  approximated  by 
an  exponential  function  of  pressure  [9,  10].  The  linear-exponential  model  equation  is  a 
piecewise  continuous  function  with  the  compliance  abruptly  changing  its  magnitude  at  the 
intersects  of  the  linear  and  the  exponential  regions.  Venegas,  Harris  and  Simon  [8],  on 
the  other  hand,  showed  that  a  single  continuous  function  in  a  form  of  sigmoidal  (tangent 
hyperbolic)  equation  represents  various  p-V  curves  extremely  well.  Parameters  in  model 
equations  (both  piecewise-continuous  and  continuous  equations)  are  determined  from  sta¬ 
tistical  processing  of  clinical  data.  More  recently  the  clinical  usefulness  of  the  sigmoidal 
model  equation  over  piecewise-continuous  representations  is  also  reported  by  the  same 
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group  of  researchers  [11]. 

Accurate  and  quantitative  determinations  of  the  form  of  p-V  model  equation  and 
its  parameters  are  prerequisite  to  the  clinical  interpretations  of  p-V  curves,  including  an 
establishment  of  ventilator  strategy  with  the  p-V  curve  guidance  in  intensive  care  for 
patients  with  acute  lung  injury  as  well  as  its  more  severe  form,  acute  respiratory  distress 
syndrome  (ARDS)  [12  -  14].  Our  objective  is  to  test  a  hypothesis  that  the  continuous 
p-V  model  equation,  particularly  in  a  form  of  an  error  function  equation,  is  effective  in 
representing  p-V  curves  from  different  sources,  and  of  different  respiratory  conditions  (a 
group  of  patients  with  ARDS  [8]  and  a  group  of  healthy  adults  [9]).  The  former  covers  both 
inflation  and  deflation  processes,  and  the  latter  includes  the  inflation  p-V  curves  before  and 
after  alveolar  recruitment  maneuver  with  a  total  of  fifty  p-V  curves.  The  report  examines 
differences  and  similarities  (1)  between  patients  with  ARDS  and  healthy  adults,  (2)  among 
patients  with  ARDS  as  well  as  among  healthy  adults,  both  in  terms  of  parameters  of  the 
error  function  p-V  model  equation. 

Equations  for  Quasi-Static  p-V  Curves 

A  model  equation,  originally  proposed  by  Venegas,  Harris  and  Simon  [8]  and  subse¬ 
quently  shown  to  represent  p-V  curves  well  for  both  inflation  and  deflation  processes  [11], 
has  the  following  sigmoidal  (tangent  hyperbolic)  form; 

dV  v  —  Vi 

—  =  -a(V-Vu)(V- VL),  -^  =  1  l+exp(-aAV(p-p0))]-'.  (ia) 

where  AV  =  Vy  -  Vi ,  Vjj  =  upper  volume  asymptote,  Vi  =  lower  volume  asymptote,  a  = 
positive  constant  and  po  =  pressure  at  the  midpoint  (inflection  point)  of  the  curve.  The 
corresponding  non-dimensional  form  of  the  sigmoidal  equation  is  [15], 

dV  A  _ 2  _  pu  —  p~w 

_  =  -_(]/  -i),  v  =  ^—z  (=  tank  (u) )  (lb) 


where 


V  =  —  ~  —  Yklll  - 

AV/2  ’  2  ’ 


1,  A  =  apo  AV. 
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Venegas,  Harris  and  Simon  suggested  that  a  p-V  equation  in  terms  of  the  error  function 
is  also  effective  in  representing  p-V  data  [8]. 

The  error  function,  erf(x ),  is  defined  as 


2  rx  2 

erf(x)  =  -j=  /  e_t  dt  with  erf{ oo)  =  1,  erf(-x)  =  -erf(x). 

V™  J o 


The  error  function  model  equation  may  be  expressed  as 


dV  _  ajAVy 
dp  4 


exp 


:-(^oAV(p  -  Po)) 


2i  YazXl 

AV 


dV  A  /  7T  2\  T7  ./a/tT  \ 
_  =  V  =  erf(—  u). 


=  i[l+er/(-AAp)],(2o) 

(26) 


Fig.l  is  a  sketch  of  a  typical  p-V  model  equation  (either  the  sigmoidal  or  the  error-function 
model  equation).  The  curve  varies  smoothly  between  the  low  pressure  asymptote,  Vl,  and 
the  high  pressure  asymptote,  Vjj.  The  midpoint  of  the  curve  where  the  volume  is  equal  to 
(Vu  +  VcJ/2  is  the  inflection  point  with  its  pressure  denoted  by  Pq.  Both  the  sigmoidal  and 
the  error-function  model  equations  are  antisymmetric  with  respect  to  the  inflection  point; 


that  is,  V(p  —  po)  —  V(po)  =  — (V(po  —  p)  —  V (po))  or  V (p)  =  —V (— p).  The  compliance, 
dV/dp ,  increases  along  the  p-V  equation  as  pressure  increases,  until’ the  inflection  point 
(=  the  point  of  maximum  compliance)  is  reached.  Then  the  compliance  decreases  with  a 
further  increase  in  pressure.  A  tangent  to  the  model  equation  curve  at  the  inflection  point 
has  the  compliance  of  c*(AV)2/4.  The  two  points  of  intersection  between  the  tangent 
and  the  two  volume  asymptotes,  V  =  Vjj  and  V  =  Vi  are  referred  to  as  the  upper  and 
lower  corner  pressure,  pCu(ci) ,  respectively.  The  pressure  difference  between  the  two  corner 
pressures  is  defined  as  the  (volume-)  gradient  pressure  range,  pgTad ■  Also,  the  pressure  at 
the  point  of  maximum  compliance  increase  (decrease)  of  the  p-V  curve,  pmci  ( pmcd ),  may 
be  specified  as  the  points  where  the  third  derivative  of  V  with  respect  to  p  is  zero. 

For  both  the  sigmoidal  and  the  error-function  model  equations, 


Pgrad 

PQ 


AV  \  _  4 
Po(dV/dp)maxJ  ~  A’ 


Pcu(cl)  ( 


Pcu(cl) 

Po 


-  V  -  <->r 


(3,4) 
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On  the  other  hand, 


_  /  Pmcd(mci)  v  _  J  (-)1.317/A  for  sigmoidal  equation; 

,nwlmu'  1  Po  \  (-)l. 596/A  for  error-function  equation.  '  ^ 

Fig.2  is  the  (non-dimensional)  p  —  V  curve,  corresponding  to  the  p-V  curve  of  Fig.l. 
The  origin  (p  =  0,  V  =  0)  represents  the  point  of  antisymmetry  (p0,  {Vu  +  VL)/2)  of  Fig.l. 
The  non-dimensional  pressure,  p,  is  the  pressure  difference,  p  -  p0,  as  a  fraction  of  p0 
(Eq.(lb).  The  normalization  of  volume  shifts  the  upper  and  the  lower  volume  asymptotes, 
Vu  and  VL  in  Fig.l  into  +1  and  - 1  respectively  in  Fig.2.  With  both  the  location  of  pQ 
and  the  volume  asymptotes  made  common  to  all  p-V  curves,  the  resulting  non-dimensional 
representations  characterize  p-V  curves  in  general  in  terms  of  a  single  non-dimensional 
parameter,  A.  (Eqs.(lb,2b))  From  Eq.(3)  the  parameter,  A,  is  four  times  the  ratio  of  the 
pressure  at  the  maximum  compliance,  p$,  to  the  volume-gradient  pressure  range,  pgrad- 
Since  the  compliance  is  maximum  at  the  origin,  the  first  quadrant  (V,  p  >  0)  in  Fig.2  is  a 
region 'of  decreasing  local  compliance  with  pressure;  while,  the  third  quadrant  (V,  p  <  0) 
is  a  region  of  increasing  local  compliance  with  pressure.  The  origin  (p  =  0,  V  =  0)  of 
dimensional  p-V  curves  is  transformed  into  (p  =  — 1,  V  {V  =  0))  on  a  p-V  curve;  hence, 
the  physiological  lower  limit  of  p  is  -1.  Various  pressure  locations  on  p-V  diagram  are 
proportional  to  1/A  as  shown  in  Eqs.(3-5).  Eqs.(3-5)  also  imply  over  the  pressure  range 
°f  P  >  0  that  there  is  no  lower  corner  pressure  (i.e  1  +pcl  <  0)  if  A  <  2,  and  that  there 
is  no  pressure  for  maximum  compliance  increase  (i.e.  1  +pmci  <  0)  if  A  <  1.317(1.596) 
for  the  sigmoidal  (error- function)  model  equation.  Both  the  sigmoidal  and  the  error- 
function  model  equations  are  capable  of  representing  p-V  curves  over  their  entire  ranges 
as  continuous  functions. 

Piecewise-continuous  model  equations  are  also  used  to  represent  p-V  curves.  Shown 
below  is  a  three-region  model  equation  [9],  relevant  to  the  present  study,  consisting  of  a 
linear  midregion  (Vj up  <  V  <  Vyip)  and  two  exponential  regions  at  high  (Vuip  <  V)  and 
low  pressure  ( V  <  VLIp)  ranges.  (The  subscripts,  L(U)IP  (lower  (upper)  inflection  point), 
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indicate  the  points  where  the  linear  midregion  equation  intersect  with  the  exponential 
equations.) 

Vmax  -  V(p)  =  (1 Vmax  ~  VuiP )  '  exp(^^^-  YEIL  (p  _  pUIp)^ 

for  VUIP  <  V 

—  Vinsuff) 

V(p) -Vlip  =  Cun(p~PLip)  for  VLip  <V  <VUIP)  (6) 

Vmin  -V(p)  =  (Vmin  -  Vup)  •  — ^  Vmm  ( p~Plip)J  for  V  <  Vlip, 

where  Ciin  =  compliance  at  the  linear  midregion,  Vmin(maxy=  volume  asymptote  of  the 
lower  (upper)  exponential  region, 

Data  Analyses 

The  two  data  sources  with 'a  total  of  50  p-V  data  sets  consists  of  (A)  ARDS  patients 
by  Harris  et  a!  (2000)  [11],  and  (B)  healthy  adults  by  Svantesson  et  al  (1998)  [9].  For  the 
data  source  A,  p-V  data  points  were  made  available  to  us  by  the  authors.  The  data  source 
B  provides  model  parameters  of  the  piecewise-continuous  model  equation,  Eq.(6),  as  well 
as  data  ranges  for  each  data  set.  Information  on  the  data  sources  relevant  to  the  present 
study  is  summarized  below. 

Data  Source  A 

21  data  sets  of  ARDS  patients  by  Harris,  Hess,  Venegas  [11],  Original  p-V  data  points 
made  available  by  the  authors,  Inflation  and  deflation  data  in  supine  position. 

Data  Source  B 

29  data  sets  of  healthy  adults  (both  male  and  female)  by  Svantesson,  Sigurdsson, 
Larsson,  Jonson  [9],  14  data  sets  before  and  15  data  sets  after  alveolar-  recruitment 
maneuver,  Inflation  data  in  supine  position 

The  parameters  of  Eq.(6)  (Vmin,  Vmax,  Vlip,  Vuip,  Vinaufj,  Cun,  Plip,  Puip)  are 
tabulated  for  all  data  sets  in  [9]. 

Data  sets  from  the  source  A  are  analyzed  by  minimizing  the  difference  between  data 
points  and  the  model  equation  (either  the  sigmoidal  or  the, error-function  model  equation) 
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through  the  application  of  the  method  of  least  squares  to  obtain  the  parameters,  A,  A  V,p0 
and  Vu  (or  VL).  To  analyze  the  data  source  B,  ten  to  twenty  five  computational  data 
points,  depending  on  the  data  range,  axe  generated  from  Eq.(6).  Then,  the  method  of 
least  squares  is  applied  to  determine  the  parameters,  A,  pQ,  AV  and  Vu  (or  VL)  of  the 
error  function  model  equation.  Parameters  of  p-V  model  equations  are  determined  for  all 
data  sets  in  Data  Source  A  (for  both  inflation  and  deflation  data)  and  B  (for  both  before- 
and  after-recruitment  maneuver).  Discussion  beyond  the  validity  test  of  the  error  function 
model  equation,  however,  is  focused  on  the  inflation  data  sets  of  the  two  data  sources. 
Results  of  data  analyses  are  summarized  in  Table.  1  (for  Data  Source  A)  and  in  Table. 2 
(for  Data  Source  B). 

Results  and  Discussion 

Fig.3  shows  a  typical  data  set  of  an  ARDS  patient  from  the  data  source  A  as  well 
as  the  sigmoidal  and  the  error  function  model  equations,  Eqs.(la)(2a),  determined  by  the 
method  of  least  squares.  The  parameters  of  the  model  equations,  (A,  pQ,  AV,  Vu),  are 
(1.470,  13.308,  3.491,  2.750)  for  the  sigmoidal  equation  and  (1.627,  13.324,  3.156,  2.584) 
for  the  error  function  equation.  Both  equations  represent  the  data  points  well  over  the 
entire  data  range.  Substantial  differences  between  the  two  equations  occur  in  high  and 
low  pressure  regions  away  from  the  data  range  as  they  approach  different  asymptotes  of 
V  =  Vu  and  Vjr,.  It  should  also  be  noted  that  there  is  no  lower  corner  pressure  for  the  data 
set  since  A  is  less  than  2  (Eq.(4)),  and  that  the  pressure  at  maximum  compliance  increase, 
Pmci,  is  very  low  at  1.385  cmH20  for  the  sigmoidal  model  equation  and  at  0.254  cmH20 
for  the  error  function  equation  (Eq.(5)). 

An  example  of  the  analysis  of  the  data  source  B  is  shown  in  Fig.4.  The  dotted  curve 
represents  the  piecewise  continuous  equation,  Eq.(6),  with 

(Knin>  Knax>  VpiP ,  Vuip,  Vinsuff,  Cun,  PLIP,  Puip) 

=  (-2230  [mL],  5870,  1513,  2884,  4125,  157  [mL/cmH20],  14.4  [cmH20],  23.1), 
reported  in  [9].  The  solid  curve  is  the  corresponding  error  function  model  equation,  Eq.(2a), 
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obtained  by  applying  the  method  of  least  squares  to  computational  data  points  generated 
over  the  data  range  of  0  <  V  <  Vinsuf/.  Since  “true”  data  points  are  unknown  in  the 
linear  region  of  Eq.(6),  we  imposed  a  constraint  that  the  inflection  point,  p0,  of  the  error 
function  equation  is  located  in  the  midregion  of  Eq.(6).  The  parameters  of  the  error 
function  equation  thus  determined  are  (A,  p0,  AV,  Vv)  =(2.0800,18.224,5.58750,  4.9114). 
Due  to  the  error  minimization  Eq.(2a)  is  nearly  identical  to  Eq.(6)  over  the  data  range 
indicated  by  the  two  triangle  marks;  which  is  valid  for  other  28  data  sets  from  Data  Source 
B.  A  continuous  change  of  the  compliance  (i.e.  non-linear  p-V  change)  in  the  region  near 
Po  has  been  described  previously  in  terms  of  the  sigmoidal  (tangent  hyperbolic)  model 
equation  [15].  It  should  also  be  mentioned  here  that,  of  twenty  nine  inflation  data  sets, 
Po  (the  inflection  point)  is  between  pup  and  puip  in  eighteen  data  sets,  equal  to  puip  in 
eight  data  sets,  equal  to  Plip  in  one  data  set,  and  po  =  puip  —  Plip  in  two  data  sets. 

Fig.5.  is  a  plot  of  the  error  function  equation,  Eq.(2a),  with  A(p/p0  -  l)/4  and 
(V  -  VL)/AV  as  x-  and  y-  axis  respectively.  All  (both  inflation  and  deflation)  data  points 
of  Data  Source  A  are  also  shown  in  the  figure,  confirming  very  good  agreements  with 
the  equation.  The  coefficient  of  determination,  R2,  is  0.999247,  which  is  comparable  in 
magnitude  to  that  of  the  sigmoidal  equation,  R2  =  0.9992  reported  in  [11],  thus  indicating 
that  both  the  sigmoidal  (tangent  hyperbolic)  p-V  equation,  Eq.(la),  and  the  error  function 
equation,  Eq.(2a),  are  very  effective  in  representing  quasi-static  p-V  curves.  Shown  in  Fig.6 
are  comparisons  between  the  sigmoidal  equation  and  the  error  function  equation  in  terms 
of  two  parameters  in  the  equations,  A  (Fig. 6(a))  and  po  (Fig. 6(b))  for  twenty  one  data  sets 
from  Data  Source  A.  Due  to  differences  in  functional  form  the  magnitude  of  A  is  slightly 
higher  for  the  error  function  equation  than  for  the  sigmoidal  equation.  On  the  other 
hand,  the  inflection  point,  po,  being  the  point  of  antisymmetry,  should  be  identical  in 
theory  for  both  model  equations.  Fig. 6(b)  confirms  it  as  the  magnitudes  of  po  determined 
by  the  method  of  least  squares  are  very  close  between  the  two  equations.  Differences  in 
the  magnitudes  of  A  and  po  between  the  two  continuous-function  model  equations  would 
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result  in  function-specific  values  for  such  quantities  as  Pcu(d),  Pmcd(mci)  of  Eqs.(4,5)  which 
characterize  p-V  curves,  indicating  the  importance  of  using  the  same  p-V  model  equation 
in  order  to  analyze  clinical  data  in  a  consistent  manner. 

Although  quantitative  comparisons  of  parameters  cannot  be  made  among  different 
p-V  model  equations  in  a  mathematically  rigorous  manner,  a  comparison  between  the 
linear-exponential  model  equation,  Eq.(6),  and  the  error  function  equation,  Eq.(2a),  is 
presented  in  Fig.7  in  a  form  of  pmci  vs  pLIp  in  Fig.7  (a)  and  pmcd  vs  pUIP  in  Fig.7  (b)  for 
Data  Source  B.  It  may  be  seen  that  the  parameters  from  both  equations  distinguish  two 
data  groups,  Before  recruitment  maneuver  and  After  recruitment  maneuver,  successfully 
and  also  that  pmcd  of  Eq.(2a)  and  puip  of  Eq.(6)  distinguish  the  two  groups  more  clearly 
than  Pmd  and  plip •  Values  of  Plip  (uip)  must  be  located  directly  from  p-V  curves;  while, 
Pmd  ( mcd )  are  automatically  generated  from  the  model  equation  once  parameters  of  the 
model  equation  are  determined. 

As  shown  in  Eqs.(2a,b),  when  p  and  V  are  made  non-dimensional  the  resulting  non- 
dimensional  p-V  equation,  Eq.(2b),  contains  A  (=the  ratio  of  (pcu  -pcl)  to  p0)  as  the  only 
parameter  representing  a  shape  of  p-V  curves  with  A/2  being  the  gradient,  dV/dp,  at 
the  origin  (where  p  =  0  i.e.  p  =  pQ.).  The  p  dV  work  associated  with  the  process  from 
the  initial  to  the  end-of-inflation  pressure  was  suggested  as  a  quantity  representing  the 
pressure  range  actually  covered  by  a  specified  p-V  curve  [15].  However,  in  analyzing  data 
sets  from  different  sources  we  found  that  the  end-of-inflation  pressure  (volume)  is  data- (or 
investigator-)  dependent,  and  may  not  be  appropriate  as  a  comprehensive  indicator  for 
data  interpretation.  Here  we  selected  two  volume  differences  to  distinguish  p-V  data  sets 
accounting  for  the  range  of  data  relative  to  the  entire  range  covered  by  the  p-V  equation; 
one  is  Vv  —  V(p  =  0)  as  a  volume  scale  indicating  the  total  available  volume  of  a  specified 
TRS,  and  the  other  is  V(p  =  20  cmH20)  —  V(p  =  0)  as  a  volume  scale  representing  volume 
range  covered  by  the  specified  TRS.  The  volume  at  p  =  20  cmH20  is  selected  arbitrarily; 
however,  in  all  the  data  sets  we  analyzed,  p-V  curves  were  obtained  beyond  p  =  20  cmH^O. 
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Fig-8  shows  A  vs  p0  in  (a)  and  [Vu  -  V(p  =  0)]  vs  [V(j>  =  20  cmH20)  -V(p  =  0)]  in  (b) 
obtained  from  the  error  function  model  equation  for  all  inflation  data  sets.  In  terms  of  A, 
a  wide  range  (1  <  A  <  6)  is  covered  by  ARDS  patients;  while,  the  range  of  A  for  healthy 
adults  is  ~  1.5  —  3.5.  The  alveolar  recruitment  maneuver  lowers  the  magnitude  of  A  of 
the  group  of  healthy  adults  as  a  whole.  A  similar  observation  may  be  made  on  pQ.  Both 
tyu  ~  V{j P  =  0)]  and  [V(p  =  20  cmH20)  —  V(p  =  0)]  in  Fig.8(b)  are  low  in  magnitude  for 
patients  with  ARDS.  On  the  other  hand,  the  recruitment  maneuver  shifts  the  location  of 
the  whole  group  to  the  right  in  Fig.8(b).  Two  data  sets,  No.20  and  5,  representing  extreme 
points  in  Fig.8  (a)  and  (b)  respectively,  clearly  show  they  are  quite  different  from  those  of 
healthy  adults,  if  Fig.8(a)  and  (b)  are  examined  together. 

According  to  Data  Source  'B  [9],  after  a  p-V  curve  before  the  recruitment  maneuver 
is  recorded,  the  lungs  are  inflated  to  an  airway  pressure  of  A0cmH2O  and  maintained  for 
15  s,  followed  by  six  pressure-controlled  breaths  (six  breaths/min.)  delivered  at  an  airway 
pressure  of  30  cmH20.  Then  a  second  large  insufflation  is  delivered  before  recording  a 
p-V  curve  after  recruitment  maneuver.  To  examine  the  p-V  curves  of  healthy  adults  as 
well  as  effects  of  the  alveolar  recruitment  maneuver  in  more  detail,  a  ratio  of  the  pressure 
at  the  inflection  point,  (p0  (before  maneuver) /p0  (after  maneuver)),  is  plotted  against 
a  ratio  of  A  (A  (before  maneuver)/ A  (after  maneuver))  in  Fig.9.  Each  data  point  is 
accompanied  by  two  numbers  indicating  the  data  set  number  and  his  or  her  age  in  the 
bracket  (unfilled  circle  for  male  and  filled  circle  for  female).  The  data  sets  of  the  younger 
may  be  seen  to  be  located  to  the  left  half  of  the  figure,  compared  to  the  older,  implying 
that,  for  the  healthy  young  adults,  A  after  the  maneuver  either  increases  slightly  or  remains 
roughly  the  same  as  A  before  the  maneuver.  In  order  to  discuss  implications  of  Fig.9 
further,  Figs.  10  and  11  show  the  error  function  p-V  equations  before-  and  after-maneuver 
along  with  the  corresponding  non-dimensional  p-V  equations  for  three  data  sets  in  Fig.  10 
from  the  region  of  A  ( before  maneuver)/ A  (after  maneuver)  <  ~  1.2,  and  for  three  data 
sets  in  Fig.ll  from  the  region  of  A  ( before  maneuver)  /  A  ( after  maneuver)  >  ~  1.2.  The 
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equations  are  plotted  over  the  measurement  range  covered  by  Data  Source  B.  Numerical 
values  of  these  data  sets  are  tabulated  below. 


Data  No.  (Age) 

1  (33) 

6  (25) 

7  (60) 

11  (55) 

13  (50) 

15  (58) 

A  (bef  ore) / A  (after) 

0.730 

0.982 

1.143 

1.390 

1.303 

1.625 

Po  (bef ore) /pQ  (after) 

1.726 

1.196 

1.840 

1.312 

2.083 

1.753 

Referring  to  Fig.  10,  the  high  p0  ratios  of  Data  1  and  7,  compared  to  Data  6,  are  results  of 
substantial  reduction  in  p0  after  the  maneuver  for  these  data  sets  as  observed  in  the  p-V 
equations.  The  triangular  marks  on  the  p  -  V  curves  in  Figs.10  and  11  indicate  locations 
of  p  =  20 cmH20;  hence,  onp^V  diagrams,  a  large  change  in  pQ  is  reflected  by  a  large 
shift  of  the  triangle  from  the  before-recruitment  location  to  the  after-recruitment  location. 
Different  degrees  of  changes  in  the  magnitude  of  A  over  the  alveolar  recruitment  maneuver 
for  the  three  data  sets  cannot  directly  be  observed  from  the  p-V  diagrams.  However,  on  the 
p-V  diagrams,  A/2  is  the  slope  of  p-V  equation  at  the  origin.(  See  Eq.(2b).)  Therefore, 
the  before-recruitment  (solid)  curve  lies  above  the  after-recruitment  (dotted)  curve  in  the 
third  quadrant  (p  <  0,  V  <  0)  for  Data  1  for  which  A  (before)  =  0.730  •  A  (after).  For 
Data  6  with  A  (before)  being  close  to  A  (after)  two  curves  are  nearly  identical.  In  Data  7 
the  after-recruitment  curve  lies  slightly  above  the  before-recruitment  curve  as  A  (be fore)  = 
1.143  •  A  (after).  The  data  sets  in  Fig.  11  all  have  the  two  ratios  well  above  unity  with  the 
high  A  ratios  resulting  in  the  after-recruitment  curves  to  lie  above  the  before-recruitment 
curves,  and  the  high  pQ  ratios  of  Data  13  and  15  being  reflected  in  the  large  shifts'  in 
triangles  between  the  two  curves  in  the  p  -  V  diagram. 

From  Eq.(2a)  the  maximum  local  compliance  (=  dV/dp  at  p  =  pQ)  may  be  expressed  as 

dV  A- AV 

-(p  =  Po)  =  __.  (7) 

Fig.12  is  a  plot  of  [dV/dp  (p  =  p0)  after  recruitment  maneuver]  vs  [dV/dp  (p  =  p0)  before 
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recruitment  maneuver]  of  Data  Source  B.  Changes  in  the  maximum  local  compliance  (= 
compliance  at  the  inflection  point)  are  small  between  before-  and  after-recruitment  data 
with  a  maximum  change  of  less  than  0.035  [. L/cmH20 ].  It  should  be  mentioned  again  that 
the  local  compliance  like  other  parameters  may  be  obtained  mathematically  in  a  continuous 
p-V  model  equation,  once  the  parameters  of  the  equation  are  determined  for  a  specified 
p-V  curve. 

Depicted  in  Fig.  13  are  p-V  curves  and  the  corresponding  p-V  curves  of  four  rep¬ 
resentative  data  sets  of  patients  with  ARDS,  drawn  over  their  measurement  ranges  The 
non-dimensional  p-V  curves  in  Fig.l3(b),  which,  we  believe,  are  more  useful  for  data 
examinations  and  interpretations,  yield  the  following  observations: 

(1.)  The  magnitude  of  A,  which  is  represented  by  the  slope  of  a  p-V  curve,  is  the  largest 
for  Data  20,  and  the  smallest  for  Data  4. 

(2.)  Since  the  origin  of  a  p  —  V  curve  is  the  location  where  the  local  compliance  is  the 
maximum  (i.e.  p  =  po),  Data  4  and  17  extend  well  into  the  region  of  decreasing  compliance, 
while,  the  pressure  range  of  Data  13  is  limited  to  the  region  of  increasing  compliance. 

(3.)  At  p  =  20  cmH20  (shown  as  a  triangle),  the  compliance  is  still  increasing  for  Data  13 
and  20,  close  to  the  maximum  for  Data  17  and  decreasing  for  Data  4. 

(4.)  The  two  volume  asymptotes,  Vv  and  VL,  are  transformed  respectively  into  V  =  1.0 
and  -1.0;  hence,  the  volume  range  of  Data  13  is  closer  to  the  lower  asymptote,  while, 
the  overall  volume  change  of  Data  4  is  small  compared  to  AV  (difference  between  the 
asymptotes). 

Furthermore  the  magnitude  of  A  in  Table  1  indicates  that  there  is  no  lower  corner 
pressure  for  Data  4.  For  each  inflation  data  set,  Table  1  lists  the  maximum  local  compli¬ 
ance.  Its  value  ranges  between  0.03  and  0.11  [ L/cmH20 ];  much  smaller  values  compared 
to  the  data  from  healthy  adults  shown  in  Fig.12.  The  maximum  local  compliance,  as 
shown  in  Eq.(7),  is  proportional  to  the  product  of  A  and  AV,  and  inversely  proportional 
to  pressure  at  the  inflection  point.  Since  both  A  and  po  are  roughly  in  the  same  order  of 
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magnitude  as  healthy  adults  (Fig.8),  the  factor  contributing  most  to  smaller  values  of  the 
maximum  local  compliance  for  patients  with  ARDS  is  A  V  as  evidenced  by  its  values  listed 
in  Table  1  and  2. 

Summary  and  Conclusions 

The  sigmoidal  (tangent  hyperbolic)  equation  is  known  to  represent  various  quasi-static 
p-V  curves  very  closely  [11].  In  the  present  study  it  has  been  shown  that  the  error  function 
model  equation  also  represents  quasi-static  p-V  curves  well  (Figs.4,5).  Major  parameters 
of  both  the  sigmoidal  (tangent  hyperbolic)  and  the  error  function  model  equations  are 
the  non-dimensional  compliance,  A,  the  maximum  local  compliance,  p0,  the  upper  (or 
lower)  volume  asymptote,  Vu  {Vu),  and  the  maximum  volume  available  for  inflation,  AV. 
Although  both  continuous  model  equations  are  antisymmetric  with  respect  to  p0,  the  non- 
dimensional  parameter,  A,  as  well  as  two  volume  asymptotes  slightly  differ  between  the 
two  equations  as  those  function-specific  parameters  are  selected  to  follow  a  specified  p- 
V  curve  as  closely  as  possible  (Fig.6).  Two  inflation  data  sources,  patients  with  ARDS 
(Data  Source  A)  and  healthy  adults  (Data  Source  B),  are  analyzed  in  detail  using  the  error 
function  p-V  model  equation  with  the  following  results; 

1.  The  alveolar  recruitment  maneuver  lowers  the  pressure  at  the  maximum  compliance, 
p0;  while,  A  remains  roughly  the  same  or  decreases  in  magnitude  (Figs.8(a),  9).  It  also 
reduces  the  upper  volume  asymptote,  Vu,  substantially  (Fig.8(b)).  The  combined  effects 
of  these  parametric  changes  due  to  the  maneuver  extend  the  range  of  p-V  curves  after  the 
maneuver  further  into  the  region  beyond  the  location  of  the  maximum  compliance  (Figs.  10, 
11). ' 

2.  The  range  of  pQ  and  A  covered  by  the  patients  with  ARDS  is  wider  than  the  corre¬ 
sponding  ranges  of  healthy  adults  (Fig.8(a)).  Substantially  lower  magnitudes  of  the  upper 
volume  asymptote,  Vu,  and  the  actual  volume  change  (Fig.8(b))  result  in  lower  values  for 
the  maximum  local  compliance  (Table  1)  compared  to  that  of  healthy  adults  (Fig.12). 

3.  The  non-dimensional  p-V  curves  combined  with  the  magnitudes  of  the  four  parame- 
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ters  of  the  model  equation  are  shown  to  help  understand  quantitatively  the  effects  of  the 
recruitment  maneuver  as  well  as  differences  among  patients  with  ARDS  (Figs.  11,  13). 

4.  An  important  advantage  of  the  continuous  model  equations  is  that  various  parameters 
characterizing  the  shape  and  the  range  of  p-V  curves,  such  as  the  maximum  local  compli¬ 
ance,  the  pressure  at  the  maximum  local  compliance  and  the  upper  and  the  lower  corner 
pressures,  may  be  evaluated  readily  from  the  model  equation  once  the  parameters  of  the 
equation  are  determined  from  p-V  curve  data.  The  mathematically  exact  relations  among 
the  parameters  also  implies  that  the  magnitude  of  either  pressure  or  volume  at  a  certain 
location  along  a  p-V  curve  may  only  be  interpreted  correctly  when  compared  to  a  char¬ 
acteristic  pressure  or  volume  of  the  p-V  curve,  as  demonstrated,  for  example,  in  Figs.  11 
and  13  when  we  discussed  the  location  of  p  -  20  cmH20  relative  to  the  pressure  at  the 
maximum  local  compliance. 
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Captions  for  Tables  and  Figures 


Table  1.  Parameters  of  Error  Function  Equation  for  Data  Source  A 

Table  2.  Parameters  of  Error  Function  Equation  for  Data  Source  B 

Fig.  1.  Continuous  p-V  model  equation. 

Fig.  2.  Non-dimensional  p-V  curve  corresponding  to  Fig.l. 

Fig.  3.  Example  of  p-V  curve  from  Data  Source  A.  filled  circle  =  Data  D, 

solid  =  error  function  p-V  equation,  dotted  =  sigmoidal  p-V  equation.  ‘ 

Fig.  4.  Example  of  p-V  curve  from  Data  Source  B.  Dotted  curve  =  piecewise  continuous 
equation,  Eq. (6),  Solid  curve  =  error  function  equation,  Eq. (2a), 

Triangle  =  upper  and  lower  data  limits. 

Fig.  5.  (V  —  Vl)/AV  vs  (^/k/4)  A(p/po  —  1)  of  Data  Source  A. 

Unfilled  .circle  =  inflation,  Filled  circle  =  deflation. 

1  - 

Fig.  6.  (a)  As  (sigmoidal  model  equation)  vs  Ae  (error  function  model  equation), 

(b)  poa  (sigmoidal,  equation)  vs  poe (error  function  equation), 
for  inflation  data  sets  of  Data  Source  A. 

Fig.  7.  (a.)  pLip[cmH20 ]  vs  pmci  [cmH^O] 

(filled  (unfilled)  circle  =  before  (after)  recruitment  maneuver). 

(b.)  puip  [cmH20]  vs  pmcd  [cmi^O]  of  Data  Source  B. 

(unfilled  (filled)  circle  =  before  (after)  recruitment  maneuver). 

(Pmci  and  Pmcd,  evaluated  from  error-function  equation.) 

Fig.  8.  (a.)  A  vs  p0,  (b .)[Vv  -  V(p  =  0)]  vs  [V(p  =  20  cmH20)  -  V{p  =0)] 
for  inflation  data  sets  from  Data  Source  A  and  B. 

Square  =  patients  with  ARDS,  Cross  =  healthy  adults  before  recruitment  maneuver, 
Triangle  =  healthy  adults  after  recruitment  maneuver. 

Fig.  9.  po  (before  maneuver)/pQ  (after  maneuver ))  vs  A  (before)/A  (after) 

Two  numbers  are  Data  Set  No.,  followed  by  his  or  her  age  in  the  bracket. 
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unfilled  circle  =  male,  filled  circle  =  female. 

Fig. 10.  p  —  V  curve  and  the  corresponding  p-V  curve  of  error  function  model  equation. 
Data  No.l  (top),  No. 6  and  No. 7  (bottom)  of  Data  Source  B  (healthy  adults). 

Solid  =  Before  recruitment  maneuver,  Dotted  =  After  recruitment  maneuver. 
Triangle  =  location  of  p  =  20  [cmi?20]. 

Fig. 11.  p  —  V  curve  and  the  corresponding  p-V  curve  of  error  function  model  equation. 
Data  No.  11  (top),  13  and  15  (bottom)  of. Data  Source  B. 

Solid  =  Before  recruitment  maneuver,  Dotted  =  After  recruitment  maneuver. 

.  Triangle  =  location  of  p  =  20  [cmH^O]. 

Fig. 12.  dV/dp  (p  =  pq)  (maximum  local  compliance)  before  recruitment  maneuver 
vs  dV/dp  {p  —  Pq)  after  recruitment  maneuver  of  Data  Source  B. 

Fig. 13.  (a)  p  —  V  curve  and  (b)  the  corresponding  p-V  curve  of  error  function  model 
equation  for  Data  No.4,  13,  17  and  20  of  Data  Source  A  (patients  with  ARDS). 
Triangle  =  location  of  p  =  20  [cxnH^O]. 
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Table  1.  Inflation  Parameters  of  Error  Function  Equation  for  Data  Source  A 


Data  No 

A 

Po 

Ay 

Vl 

(dV/dp)  at  po 

1. 

2.9652, 

22.39, 

2.3726, 

-0.0709 

0.0785 

2. 

2.7173, 

21.98, 

1.5612, 

-0.0762 

0.0482 

3. 

3.3532, 

25.08, 

1.6193, 

-0.0365 

0.0541 

4. 

1.6273, 

13.32, 

3.1567, 

-0.5727 

0.0963 

5. 

2.7160, 

30.36, 

1.6216, 

-0.0768 

0.0362 

6. 

2.6029, 

23.88, 

1.5066, 

-0.0962 

0.0410 

7. 

1.7288, 

20.15, 

3.0989, 

--0.3645 

0.0664 

8. 

1.8901, 

14.95, 

1.7905, 

-0.2185 

0.0565 

9. 

3.5379, 

25.24, 

2.7570, 

-0.0350 

0.0965 

10. 

2.7981, 

26.20, 

3.7326, 

-0.1583 

0.0996 

11. 

S 

2.4296, 

17.89, 

1.3424, 

-0.0997 

0.0455 

12. 

1.2500, 

11.59, 

1.2470, 

-0.2829 

0.0336 

13. 

3.2725, 

29.86, 

3.9075, 

-0.0861 

0.1070 

14. 

2.7915, 

15.31, 

1.6284, 

-0.0705 

0.0742 

15. 

1.9487, 

18.37, 

1.8797, 

-0.2423 

0.0498 

16. 

2.4566, 

26.98, 

1.3277, 

-0.0787 

0.0302 

17. 

2.1027, 

19.31, 

1.3316, 

-0.1279 

0.0362 

18. 

1.1672, 

13.98, 

2.0685, 

-0.5149 

0.0431 

19. 

3.1381, 

26.80, 

3.1306, 

-0.0851 

0.0916 

20. 

5.4709, 

30.03, 

1.7695, 

-0.0097 

0.0805 

21. 

3,2818, 

24.43, 

2.8956, 

-0.2075 

0.0972 

Po  in 

[cmi720],  AV"  and  Vj_,  in  [L],  dV/dp  in  [L/cmi720]. 
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Table  2.  Parameters  of  Error  Function  Equation  for  Data  Source  B 


Data  No 

A 

Po 

AV 

vL 

Vv 

l.F 

B. 

1.5136 

25.00 

5.1043 

-0.8889 

4.2154 

A. 

2.0711, 

14.48 

2.7130 

-0.2700 

2.4430 

2.r 

B. 

2.5520 

31:70 

5.4303 

-0.2863  . 

5.7166 

A. 

3.F 

1.7587 

18.17 

4.1416 

-0.6590 

3.4826 

B. 

2.0979 

21.76 

3.7676 

-0.3399 

3.4277 

A. 

2.2960 

18.19 

3.1331 

-0.2300 

2.9031 

4.M 

B. 

2.1491 

21.11 

4.5443 

-0.3726 

4.1717 

A. 

5.F 

2.1887 

18.57 

4.4023 

-0.3554 

4.0469 

B. 

2.2192 

24.60 

6.4953 

-0.4971 

5.9982 

A. 

6.M 

1.5729 

13.46 

5.2884 

-0.8895 

4.3989 

B. 

2.3950 

23.21 

4.2092 

-0.2567 

3.9525 

A. 

7.M 

2.4387 

19.40 

3.9049 

-0.2722 

3.6327 

B. 

2.6997 

31.00 

7.3180 

-0.2167 

7.1013 

A. 

8.F  1 

2.3615 

16.85 

5.5288 

-0.5641 

■  4.9647 

a 

2.2901 

30.80 

8.2489 

-0.5185  . 

7.7304 

A. 

9.F 

1.9974 

17.20 

5.0638 

-0.5478 

4.5160 

B. 

— 

— 

- — 

•  —  — .  _ 

A. 

10.M 

2.0340 

15.90 

3.7392 

-0.4068 

3.3324 

B. 

3.4081 

28.80 

7.211- 

-0.0491 

7.1627 

A. 

ll.M 

2.3998 

17.92 

5.7276 

-0.5502 

5.1774 

B. 

3.1746 

25.86 

7.0726 

-0.1361 

6.9365 

A. 

12.F 

2.2828 

19.70 

6.79820 

-0.8620 

5.9362 

B. 

2.5060 

23.20 

5.0810 

-0.2476 

4.8334 

A. 

13.F 

2.3844 

16.36 

4.3799 

-0.3799 

4.0000 

B. 

2.3600 

29.17 

6.81755 

-0.4353 

6.3822 

A. 

14.M 

1.8107 

14.00 

4.03567 

-0.6883 

3.3472 

B. 

2.2260 

24.55 

6.7170 

-0.4828  . 

6.2342 

A. 

15. F 

2.0800 

18.22 

5.5875 

-0.6761 

4.9114 

B. 

3.5047 

25.70 

4.5022 

-0.0027 

4.4995 

A. 

2.1561 

14.65 

3.3016 

-0.4836 

2.8180 

B  —  Before  recruitment  maneuver,  A  =  After  recruitment  maneuver. 
Po  in  [cmH20],  AV,  VL  and  Vv  in  [L]. 
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Fig. 3.  Example  of  P-V  curve  from  Data  source  A. 
filled  circle  =  Data  D,  solid  =  error  function  p-V  equation,  dotted  =  sigmoidal  p-V 
equation 
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Fig. 4.  Example  of  p-V  curve  from  Data  Source  B. 

Dotted  curve  =  piecewise  continuous  equation,  Eq(6),  Solid  curve  =  error  function  equation 
Eq(2a),  Triangle  =  upper  and  lower  data  limits. 
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Fig.6.  (a)As  (sigmoidal  model  equation)  vs  Ae  (error  function) 
model  equation. (b)  p0  s  (sigmoidal  equation)  vs  pQ  e  (error  function 
equation)  for  inflation  data  sets  of  Data  Source  A. 
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Pmcd  Pmci 


Fig.7.  (a)  pup  [cmH20]  vs  pmc,  [cmH20](filled  (unfilled)  circle=  before 
(after)  recruitment  maneuver),  (b)  py//>[cmH20]  vs  pmcd  [cmH20]  of  Data 
Source  B. (unfilled  (filled)  circle  =  before  (after)  recruitment  maneuver.  pmcj 
and  Pmcd  evaluated  from  errorfunction  equation.) 
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6 


Po  [cm  H2O] 


V(P=20)-V(P=0)  [L] 

Fig.8.  (a.)  A  vs  p0  ,  (b.)  \yu  -V(p  =  0)]  vs  \V(P  =  20cmH2O)-V(P  =  0)] 
for  inflation  data  sets  from  Data  Source  A  and  B. 

Square  =  patients  with  ARDS,  Cross  =  healthy  adults  before  recruitment  maneuver 
Triangle  =  healthy  adults  after  recruitment  maneuver. 
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Fig.  9.  p0  (before  maneuver/ p0 after  maneuver  vs  A  (before)/ A  (after). 

Two  numbers  are  Data  Set  No.,  followed  by  his  or  her  age  in  the  bracket, 
unfilled  circle  =  male,  filled  circle  =  female. 


jaAnaueifl  J9UV  *xb|/\|  (dp/Ap) 
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Chapter  2.  A  Mechanistic  Model  for  Quasi-Static  Pulmonary  Pressure- Volume 
Curves:  Model  Development  for  Inflation  Process 

Abstract 

A  mechanistic  model  of  a  total  respiratory  system  is  proposed  to  understand  differences 
in  quasi-static  pressure- volume  (p  -  V )  curves  of  the  inflation  process  in  terms  of  the  alve¬ 
olar  recruitment  and  the  elastic  distension  of  the  wall  tissues.  In  the  model,  based  on  the 
Boltzmann  statistics,  the  total  respiratory  system  consists  of  a  large  number  of  elements, 
each  of  which  is  a  subsystem  of  a  cylindrical  chamber  fitted  with  a  piston  attached  to 
a  spring.  The  alveolar  recruitment  is  simulated  by  allowing  a  distribution  of  the  critical 
pressure  at  which  an  element  opens;  while  the  wall  distension  is  represented  by  the  piston 
displacement.  Various  parameters  in  the  error-function  p  —  V  model  equation  are  related 
to  the  properties  of  the  mechanistic  model.  The  parameters  of  the  model-based  p-V  equa¬ 
tion  are  determined  for  each  clinical  data  set  for  a  total  of  twenty  one  p  —  V  data  sets  of 
patients  with  acute  respiratory  distress  syndrome  by  a  computational  minimization  proce¬ 
dure  between  the  equation  and  the  data  points,  results  of  which  show  excellent  agreements 
between  the  two. 
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Nomenclature 


A 

B 

C 

f,  F 

Ii  (*  =  1  -  5) 
k 
N 
N, 

P 

P 

Pcj 

Pcj 


{k/As)y0  =  p0 

piston  surface  area  on  which  pressure  is  acting  (Fig.4) 

(k/A3)yT  =  pQ  •  yT o,  threshold  pressure  for  onset  of  saturation 

(P/2)1/2  -pQ  =  V^A/4 

distribution  functions  (Eqs.(6)  k,  (13c)) 

functions  defined  in  Eq.(13b) 

spring  constant  [N/m]  (Fig.4) 

total  number  of  TRS  elements 

number  of  elements  at  energy  level  j 

pressure  (interpleural  pressure  difference) 

non-dimensional  pressure,  p/po  —  1 

critical  pressure  at  which  an  element,  j,  ‘pops  open’. 

Pcj/Po  (Eq.(13c)) 


PS 

Po 

U(p) 

AU 

V 

VP 

VL{U) 

V 

Vj 


pressure  at  the  end  of  inflation  process 

pressure  at  the  inflection  point  in  p-V  equation,  Eq.(l) 

total  energy  of  TRS  at  p=p 

=  U(p  =  pf )  -  U(p  =  Q) 

volume 

volume  change  from  the  state  of  p  =  0 
lower  (upper)  volume  asymptote  (Fig.l) 
non-dimensional  volume,  (V  -  ((Vu  +  VL)/2)/(AV/2),  (Eq.(2)) 
volume  of  an  element  j 


Ay  Vu  -  VL  =  NV0(yT0  +  l) 

Vo  ‘pop-open’  volume  (=  A3y0)  (Fig.4) 

Vj  piston  displacement  of  an  element  j  (Fig.4) 
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2/o  ‘pop-open’  displacement,  (=  Vq/As)  (Fig.4) 

Vt  piston  stroke  limit  (Fig.4) 

Vto  Vt/Vo 

Greek  symbols: 

constant  of  proportionality  (Eq.(la)) 
constant  in  the  distribution  function  (Eq.(4)) 
=  {A\/k)0{  Eq.(5)) 
energy  stored  in  an  element  j  (Eq.(3)) 
apo AV(non-dimensional  parameter)  (Eq.(l)) 
(8/7t)  1/2/ A,  Standard  deviation  (Eq.(13c)) 
Ap/2  (Eq.(2))  . 

Superscript: 

related  to  a  single  TRS  element 

Acronyms: 

ARDS  acute  respiratory  distress  syndrome 

TRS  total  respiratory  system 


a 

0 

P 

A 

a 

u 
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Introduction 


Quasi-static  pulmonary  p-V  (pressure  -  volume)  curves  are  used  routinely  to  obtain 
quantitative  information  on  the  respiratory  system  that  is  important  for  both  research  and 
clinical  guidances,  as  the  conditions  of  gas  exchange,  the  primary  role  of  the  respiratory 
system,  are  related  to  the  characteristics  of  the  curve.  During  the  inflation  (inspiration)  and 
the  deflation  (exhalation)  processes,  the  respiratory  system  changes  its  volume  (measured 
m  L  (=  10~3  m3)  or  mL),  lung  (alveolar)  pressure  as  well  as  the  pleural  pressure  (the 
pressure  of  the  thin  liquid  film  that  couples  the  lungs  and  the  chest  wall  pleurae).  The 
pressure,  p,  refers  to  the  interpleural  pressure  difference  (i.e.  the  difference  between  the  lung 
pressure  and  the  pleural  pressure)  measured  in  water  head  [cm-H20]  (1  cmH20  =  98  Pa). 
Clinical  p-V  curves  are  commonly  obtained  for  an  anesthetized  human  subject  in  supine 
position  by  sequentially  adding  (or  withdrawing)  incremental  gas  volumes  (~  50-100  mL) 
in  a  stepwise  manner  (with  a  duration  of  ~  5  seconds  per  step)  [1,2].  Fig.l  is  a  typical 
inflation  p-V  curve,  consisting  o'f  a  nearly  linear  region  of  high  compliance  (i.e.  large 
dV/dp)  sandwiched  between  two  segments  with  low  compliance  at  low  and  high  pressure 
regions.  The  shape  of  the  p  -  V  curve  is  affected  by  two  mechanisms,  the  distension  of  the 
elastic  respiratory  wall  tissue  components  and  the  recruitment  (‘pop-open’  mechanism) 
of  the  alveoli.-  The  latter  is  the  opening  of  alveoli  overcoming  the  surface  tension  at 
the  interface  between  the  gas  and  the  liquid  film  lining  the  alveolar  surface.  A  pressure 
increase  (i.e.  an  increase  in  the  interpleural  pressure  difference)  results  in  the  recruitment 
of  a  greater  number  of  alveoli.  The  high  compliance  is  believed  to  be  associated  with  both 
the  distension  of  open  parts  of  the  lungs  and  the  (alveolar)  recruitment  of  collapsed  parts 
of  the  lungs  [3].  Some  protective  ventilation  strategies,  based  on  patients’  quasi-static  p-V 
curves,  have  been  proposed  for  lung  disease  patients  in  intensive  care  units.  Amato  and 
coworkers  [4,5]  demonstrated,  based  on  their  clinical  study  involving  patients  with  acute 
respiratory  distress  syndrome  (ARDS),  that  a  ventilator  strategy  guided  by  the  p-V  curve 
resulted  in  reduced  lung  trauma,  a  high  weaning  rate  and  improved  survival  compared  with 
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a  conventional  ventilator  strategy  without  the  p-V  curve  guidance.  Also,  a  recent  ARDS 
Network  report  [6]  on  a  clinical  study  involving  861  patients  shows  lower  mortality  in  the 
group  treated  with  lower  tidal  volume  than  in  the  group  treated  with  traditional  higher 
tidal  volumes.  Although  a  use  of  p-V  curves  is  not  mentioned,  the  report  underscores  the 
importance  of  optimized  ventilator  strategy. 

In  order  to  quantify  the  characteristics  of  p-V  curves  as  well  as  their  changes  observed 
in  clinical  settings,  various  p-V  model  equations  have  been  proposed  [7-12].  Parameters  in 
model  equations  are  determined  from  statistical  processing  of  clinical  data.  It  is  important 
that  these  parameters  should  have  some  physiological  interpretations.  Also,  to  understand 
the  shape  of  p-V  curves  in  terms  of  mechanical  behavior  of  lungs,  multi-compartment  lung 
models  were  developed  and  used  to  obtain  information  on  the  effects  of  lung  elasticity  and 
a  degree  of  alveolar  recruitment  on  p-V  curves  [3,13].  Although  these  analyses  serve  to 
relate  the  internal  elastic  conditions  of  the  total  respiratory  system  (TRS)  to  general  p-V 
curve  behavior,  there  has  been  no  .attempt  to  interpret  individual  differences  in  p-V  curves 
directly  in  terms  of  internal  elastic  properties,  alveolar  recruitment  and  their  changes. 
Prom  an  analytical  viewpoint  the  quasi-static  p-V  curves  are  more  amenable  to  theoretical 
studies  because  at  each  state  we  may  be  able  to  apply  equilibrium  principles.  An  overall 
objective  of  this  report  is  to  test  the  hypothesis  that  a  mechanistic  model,  based  on  the 
continuous  alveolar  recruitment  and  the  elastic  distension  of  the  wall  tissues,  is  effective  in 
understanding  a  relation  between  the  observed  pulmonary  behavior  (as  p-V  curves)  and 
the  corresponding  internal  respiratory  response  (in  terms  of  the  mechanistic  model). 
Continuous  Equation  for  Quasi-Static  p-V  Curves 

In  the  past  piecewise-continuous  equations  were  used  to  generate  such  quantities  as 
compliance,  the  lower  and  upper  inflection  points  that  may  reflect  the  internal  conditions 
of  TRS  [7-12].  There  are  two  continuous  model  equations  that  simulate  various  p-V  curves 
accurately  over  the  entire  range  of  p-V  data.  One  is  a  sigmoidal  (tangent  hyperbolic) 
equation,  and  the  other  an  error-function  equation,  both  originally  proposed  by  Venegas 
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and.  his  coworkers  [1].  Since  the  analytical  development  to  follow  utilizes  the  error  function 
representation  of  p-V  curves,  its  characteristics  are  discussed  below  in  some  detail. 

The  error  function  p-V  equation,  plotted  in  Fig.l,  may  be  expressed  as 

»)'(£->)'].  o.) 

■"(£»(£-));  m 

where  AV  =  V\j  —  Vl,  Vjj  =  the  upper  asymptote,  Vl  =  the  lower  asymptote,  a  = 
positive  constant,  po  =  a  pressure  at  the  midpoint  (inflection  point)  of  the  curve  and  A 
(non-dimensional)  =  ap0  AV  [1,2,14]  with 


e  t2dt  erf( oo)  =  1,  erf(-x)  =  -erf(x). 

The  corresponding  non-dimensional  forms  are, 

dV  A  ,  7T  ox  77  wV7T  N  x 

,  ^  =  2  •^’(-4“)-  V  =  er/(^-“)-  (2) 

where  V  =  .[V  -  (Vv  +  VL)/2]/(AV/2),  p  =  (p/p0)  -  1,  w  =  Ap/2. 

Eq.(2)  satisfies  the  following  conditions:  „ 

V  (p  =  0)  =0,  V  (p)  =  —V  (— p)  (antisymmetry  with  respect  to  p  =  po), 

dV/dp  (p  =  0)  =  A/2,  F  (p  ±00)  =  ±  1,  dV  (p  -4  ±oo)/dp  =  0. 

A  clinical  data  source  of  p-V  curves  we  use  in  the  present  analyses  are  twenty  one  data 
sets  of  ARDS  patients  (both  inflation  and  deflation  data)  in  supine  position  by  Harris  et 
al  [2],  made  available'  by  the  authors.  Data  sets  are  analyzed  by  minimising  the  difference 
between  data  points  and  the  error  function  model  equation  through  the  application  of  the 
method  of  least  squares  to  obtain  the  parameters,  A,  AV,p0  and  Vy  (or  VL).  Plotted  in 
Fig. 2  (Fig. 3)  are  264  inflation  data  points  (225  deflation  data  points)  in  terms  of  Eq.(lb), 
(V  —  Vl)/ AV  vs  v/7rA(p/po  -  l)/4.  Agreement  is  excellent  between  the  data  and  the  p-V 
equation  with  B?  (the  coefficient  of  determination)  =  0.99938  for  the  inflation-  and  = 
0.99907  for  the  deflation-  data  points. 
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Development  of  a  Mechanistic  Model 

An  overall  objective  of  the  development  of  a  mechanistic  TRS  model  is  to  derive  an 
equation  for  the  volume  participating  in  the  p-V  variations.  We  consider  a  TRS  comprised 
of  a  very  large  number  of  elements  with  N  =  total  number  of  elements.  Based  on  the  char¬ 
acteristics  of  a  single  element  that  are  common  to  all  elements,  a  distribution  of  elements 
is  derived  over  a  distribution  parameter.  The  mechanistic  model  of  an  element  is  shown 
in  Fig.4.  An  arbitrary  element,  j,  consists  of  a  cylindrical  chamber  containing  a  piston 
(with  its  surface  area,  Aa,  [m2]),  which  is  attached  to  a  spring  with  its  spring  constant, 
k  [N / m] .  The  element  is  closed  when  the  piston  is  located  at  the  left  end  of  the  cylinder  in 
Fig.4.  When  pressure  acting  on  the  left  end  of  the  piston  reaches  a  certain  critical  value, 
the  piston  suddenly  moves  to  a" new  position  (‘pop-open’  mechanism)  with  the  elemental 
volume,  Vo,  in  the  figure  indicating  an  elemental  volume  increase  due  to  the  sudden  piston 
displacement  of  y0 •  Once  the  element  is  open  with  its  volume  of  Vo,  any  further  increase  in 
pressure  results  in  a  volume  increase  as  the  piston  moves  to  the  right  until  it  reaches  the 
end  of  the  cylinder.  (The  symbol, indicates  an  elemental  quantity.)  In  the  model  the 
pop-open  volume,  Vo  (=  Aayo),  and  a  further  volume  increase  due  to  piston  displacement 
represent  the  opening  of  alveoli  and  the  elastic  distension  of  the  wall  tissues  respectively. 
We  define  pc j  as  the  critical  pressure  at  which  the  element,  j,  ‘pops  open’. 

Referring  to  Fig.4,  the  elemental  volume,  V),  at  p  (>  pcj)  is  equal  to  %  +  Aayy,  which, 
upon  application  of  a  quasi-static  force  balance  across  the  piston,  Aa  ( p—pCj )  =  ky-j,  may 
be  expressed  as,  Vj  =  Vq  +  {A2a/k)  (p—pCj)-  Also  the  piston  position  of  an  element  reaches 
its  stroke  limit  of  yx  when  pressure,  p,  reaches  (pc  j  +  {k  j  Aa)yx)  ■  The  mechanistic  model  of 
an  element,  therefore,  goes  through  three  stages  in  the  inflation  process  —  closed,  open  & 
unsaturated  (i.e.  y3  <  yT)  and  open  &  saturated  (i.e.  y3  =  yT)-  The  model  assumes  that 
the  critical  ‘pop-open’  pressure,  pCj,  as  well  as  the  location  of  the  piston  for  open  elements, 
Vj,  vary  from  element  to  element  at  an  arbitrary  quasi-static  state  (p,  V),  and  that  other 
quantities  such  a s  k,  Aa,  Vq  and  yx  are  constant  and  common  for  all  elements.  The  energy 
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level  of  an  open  and  unsaturated  element,  j,  consists  of  the  activation  energy  required  to 
pop  open  the  element,  tj  A,  and  the  energy  stored  in  the  spring,  tjS\  For  tj  a,  we  assign 
the  compression/ expansion  work  under  constant  pressure;  i.e.  tjA  =  pCjV{ j;  while,  tj  5  is 
equal  to  ky‘j/2,  which  may  be  expressed  in  terms  of  pcj  as  tjs  =  (A23/2k)(p  -  pCj)2  from 
an  application  of  the  force  balance. 

In  summary,  at  a  quasi-static  (-equilibrium)  state  at  p  =  p,  a  TRS  element,  j,  belongs  to 
one  of  the  following  states: 

if  P  <  Pc  j,  the  element,  j,  is  closed  with  Vj  =  0,  tj  =  0. 

if  pCj  —  p,  the  element,  j,  pops  open  with  Vj  =  Vq,  tj  =  pCjVo. 

if  p  —  ( k/As)yT  <  pCj  <  p,  the  piston  of  the  open  and  unsaturated  element,  j, 
moves  to  a  location,  yjt  with  Vj  =  %  +  Aa  yjt  tj  =  pcjV0  +  (A23/2k)(p  -  pcj)2. 
if  Pcj  <  P  —  {k/As)yT,  the  piston  of  the  open  and  saturated  element,  j, 
remains  at  the  stroke  limit,  yT,  with  Vj  =  V0  +  A3yT,  tj  =  pcjV0  +  {k/2)y\. 

The  state  of  an  element  follows  the  sequence  above  during  the  inflation  process  as  p  in¬ 
creases.  To  obtain  an  explicit  form  of  the  distribution  function  of  TRS  elements  over  the 
distribution  parameter,  pcj,  we  focus  on  open  and  unsaturated  elements,  for  which  the 
elemental  energy  shown  above  may  be  rewritten  as, 

tj  {open,  unsaturated )  =  ^  \ipcj  -  (p  -  ^)] 2  +  y  (2 p-  .  (3) 

According  to  the  Boltzmann  statistical  model  ([15,16]  for  example),  which  assumes  that 
there  is  no  limit  in  the  number  of  elements  per  energy  state,  the  most  probable  distribution 
Nj/N  (a  fraction  of  elements  at  an  energy  level,  tj),  may  be  expressed  as 

constant )  (4) 

A  substitution  of  Eq.(3)  into  Eq.(4)  with  the  summation  replaced  by  an  integral  over  the 
whole  range  of  pc  j  for  a  large  number  of  elements,  yields 

J^-  =  eXp(-^(p-Pcj-±V0)^/ jjxp(-i(p-Pcj.^voyypct  (5) 


Nj/N  =  e-^ 


{(5  =  unspecifie 
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where  0  (As/k)0.  It  should  be  noted  that  the  integration  in  the  denominator  ranges 

from  -co  and  +oo.  As  summarized  above,  any  arbitrary  element  remains  active  (open 
and  unsaturated)  only  in  a  certain  range  of  pcj ;  therefore,  the  application  of  Eq.(5)  over 
the  entire  range  of  -  oo  <  pcj  <  oo  assumes  that  the  distribution  function  that  is  valid  for 
active  elements  is  also  applicable  in  evaluating  the  number  of  closed  as  well  as  saturated 
elements.  Then,  upon  performing  the  integration  of  the  denominator  in  Eq.(5),  we  obtain 
for  dNj/N(=  a  number  fraction  of  elements,  for  which  the  magnitude  ofpcj  ranges  between 
Pcj  and  pcj  +dpcj), 

=  f-  dV*ji  f  =  [7^  j  *«*(•-  §•(/>  -  Pcj  -  JjVo)2)  (6) 

Noting  that  elements,  j,  with  pCj  in  the  range  of  0  <  pc  j  <  p,  are  open  at  p  —  p,  and 
that  (kJAs)  yT  (=  B)  is  the  pressure  at  which  an  element  j  with  pcj  =  0  reaches  the  piston 
stroke  limit  of  yT,  the  volume  change  with  pressure  needs  to  be  evaluated  for  the  following 
two  pressure  ranges;  pressure  range  1:  0  <  p  <  B  and  pressure  range  2:  B  <  p. 

Pressure  Range  1:  0  <  p  <  B 

Since  the  pressure  is  below  B  (the  threshold  pressure  for  the  onset  of  saturation),  all 
open  elements  are  active  (unsaturated)  with  yj  <  yrp.  Then  the  total  volume,  Vp(=  V(p  — 

P)  ~  V(p  =  0) ;  i.e.  a  volume  change  from  the  state  of  p  =  0),  participating  in  the  inflation 
process  is, 


Vp  (p)  = 


N\  f  (Vo  +  A,yj)f(p  =  p)dpcj  +  A,yj  [  f(p  =  p)d.pcj  . 

[‘'Pcj— 0  J  Pc  j  — P 


The  first  term  on  the  right  hand  side  represents  a  volume  increase  due  to  the  elements  that 
pop-open,  followed  by  a  piston  displacement  (yj  =  A](p  -  pcj)/k  from  the  force  balance) 
during  the  inflation  process  from  p  =  0  to  p  =  p ;  while,  the  second  term,  noting 


/°°  /*CO 

f(p  —  p)  dpcj  —  /  f(p  =  0) 

rC  j~P  Jpcj=0 


accounts  for  the  elements  that  are  already  open  at  p  =0  and  the  piston  displacement 
Vj  (=  A]p/k)  is  the  only  mechanism  available  for  the  volume  increase. 
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After  expressing  y.j  in  terms  of  p  and  pc j  as  shown  above,  the  equation  may  be  written  as, 


=  iV  f- 


1  ~  A2  fl  r°°  4  2  a 

( 2yo  +  -f  t)  •  exP{-  f  *2)  dt  +  J  Yp  '  exp(~  2z^dz  ^ 


where  A  =  (k/Aa)  y0 ,  t  =  -z  =  (p- Pcj  -  k V0/A2). 

Pressure  Range  2:  B  <p 

In  this  pressure  range,  the  elements  with  0  <  pc  j  <  p  -  B  as  well  as  the  elements  that 
are  already  open  at  p  =0  are  saturated  (i.e.  yj  =  yT  for  the  elements);  while,  the  elements 
with  p-  B  <pcj  <p  remain  unsaturated  (i.e.  yj  <  yT);  therefore,  ' 


Vp  =N 


fv  B  ~  ^  rp^ 

(V0  +  AsyT)  •  f(p  =  p)  dpCj  +  N  (V0  +  Aayj)  -  f(p  =  p)  dpcj 

JQ  Jv-B 


+  N  AsyT  •  f(p  =  p)  dpCj. 


Eqs.(7)(8),  after  integration,  become, 

I  „ 

Vi  (o  <  V  <  B)  =  m  [/,  +  i  £(1  -  /,)  -  m  +  4(p)] , 

Vr(  B  <  p  )  =  W0[/1  +  ^»(i-/l)  +  ^Lzi/4__|_  +  fc2+i/: 


where 


h  =  erf(C),  I2(p)  =  exp(-C2(^  -  l)2)  -  exp(-C2),  I3(p)  =  erf(C( |  -  1)), 

U  =  erf(C(l  -  j/?o)),  /5  =  exp(-C2(yT0  -  l)2)  -  exp(-  C2), 

C  =  (ft/ 2)  V2  •  A,  yT0  =  i/T/yo  =  B/A . 

Therefore,  the  mechanistic  model  yields  the  following  p-V  equation; 


V  (0  <  p  <  B)  =  Vu 


+  ^  (^  (*  -  fa,)  -  &Lti  -  .  /4  +  hzjgt  +  /3(p))i  (10a) 


l/(B<p  )  -  yu  yVoWro  +  I)  |  NWyro  ±  1) 


h(p)- 


(106) 
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In  order  to  relate  the  present  TRS  model,  intended  to  describe  the  internal  respira¬ 
tory  conditions,  to  the  p-V  curve  that  quantitatively  describes  overall  variations  in  TRS 
conditions,  an  additional  condition  that  needs  to  be  satisfied  is  conservation  of  energy.  For 
a  quasi-static  process  from  an  initial  state  of  (p  =  0,  V  =  V(p  =  0))  to  a  final  state  of 
(P  =  Pf,V  =  V(p  =  pf))  (p f  =  the  pressure  at  the  end  of  a  measured  p-V  curve),  the 
conservation  of  energy  neglecting  any  dissipative  mechanisms  may  be  written  as, 


rV(p=pf) 

AU(=U(p  =  pf)-U(p  =  0))  =  pdV,  (11) 

J  V  (p=0) 

where  U[p)  represents  the  total  energy  of  TRS  at  p  =  p  to  be  evaluated  from  our  mecha¬ 
nistic  model;  while  the  right  hand  side  of  the  equation  is  work  associated  with  the  inflation 
process  that  must  be  evaluated  from  the  p-V  model  equation.  A  further  development  of 
the  energy  equation  will  be  discussed  in  the  next  section. 

Mechanistic  Model  vs  Error  Function  p-V  Equation 

Relations  between  the  parameters  in  the  error  function  p-V  equation  and  the  pa¬ 
rameters  in  our  mechanistic  model  are  derived  based  on  the  observation  that  the  p  -  V 
relation,  Eq.(lOb),  for  the  high  pressure  region  as  well  as  the  corresponding  equation  for 
the  local  compliance, 


v  ( B  <  P )  =  V„  -  +  V  +  EMVTO  +  D  .  er/(c  (£_!)) , 

become  identical  to  the  error  function  model  equation,  Eq.(l), 


V  =  Vu-^+(^).er/^A(^-  l)), 


dV  a 
dp 


if  we  set 


NVQ{yT0  +  1)  =  AV,  A^=  ( k/As )  jfo)  =  p0,  C  (=  1  •  a)  =  A. 


(12) 
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Before  further  developments  are  made  on  the  mechanistic  model,  our  results  are  summa¬ 
rized  below,  based  on  the  parametric  relations,  Eq.(12),  between  the  error  function  p-V 
equation  and  the  model. 

Pressure  -  Volume  (p  -  Vp)  Equation: 


Vp  (0  5:  P  <  Po  •  Pro)  = 
Vp  (Po  ■  VTO  <  P )  = 


AV 


Vto  + 
AV 


r  _L  yT0  a  ,  VTO -l  T 

VTo  +  l\h  +  ~{1~h)  +  ~2-Ia 


2  h  ,  Pro  +  1  T  , 

^A+~T~'h{ 


Pressure  -  Volume  (p  —  V)  Equation: 


V  (0  <  p  <  po  •  yro )  = 


Vu  +  VL 


+ 


AV 


Vto  + 

V ( po  •  pro  <  p) 


VU  +  VL  ,  AV 
— —  +  -J-h(p) 


where  Ix  =  erf(C),  I2{p)  =  exp{-C 2  p2)  -  exp(-C2), 

h(p)  =  erfiCp),  h  =  erf(C(  1  -  p^)), 
h  =  exp(-C2  ( pro  -  l)2)  -  exp{-  C2),  C  =  a/7tA/4. 


Distribution  Function: 


(13a) 


T  ilzii)  g, + 1  _  M  -yj^lu+ , 


(136) 


dN, 

N  ■  dpCj 


m, 


where  pcj  =  pcj/p0,  a  =  (  8/tt^/A. 


(13  c) 


The  model-based  p-V  equation,  Eq.(13b),  consists  of  two  regions.  The  solution  for  the 
high  pressure  region  is  identified  with  the  error-function  p-V  equation.  The  p-V  equation 
for  the  lower  pressure  region  as  well  as  the  boundary  pressure  between  the  high-  and 
low-pressure  region  contain  the  parameters  of  the  p-V  equation,  A,  p0,  AV,  Vy  (orVL)), 
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and  an  additional  parameter,  y? 0.  Conservation  of  energy,  Eq.(ll),  is  utilized  to  find  the 
magnitude  of  j/to ■  Similar  to  the  p-V  model  equation,  Eq.(13b),  the  evaluation  of  Eq.(ll) 
depends  on  the  magnitude  of  the  final  pressure,  pf,  relative  to  the  boundary  pressure, 
PoVro,  between  the  high  and  the  low  pressure  regions  of  the  model-based  p-V  equation. 
As  will  be  shown  later  in  the  analyses  of  clinical  data,  the  magnitude  of  yTQ  is  less  than 
unity;  hence,  the  conservation  of  energy  is  further  developed  for  the  case  of  p0-yT0  <  pf  (i.e. 
yro-l  <  P/(=  Pf/po  -  1)).  Accordingly  the  left  hand  side  of  Eq.(ll)  may  be  evaluated 
from  the  elemental  distribution  function,  F(p)  of  Eq.(13c),  along  with  the  elemental  energy 
summarized  in  the  paragraphs  preceding  Eq.(3),  yielding, 


N  Pn  A2  /■Pz  +  l-yro  i 

'  ~k  I  +  ffio) -  F(P  =  Pf)  dPcj 


rPj+l  x 

+  /_  -  (Pel  +  9  (Pf  +  1  -  Pc  j?  )  •  F(p  =  Pj)  dp Cj 

J  P  f  +  l—yTo  " 

NPoAl  [  _  _ 

l  /  KPci  -  Pf)  F(p  =  Pf)  dpCj 


Vth  fPf+l~yT0  rPf+ 1 

+  —  yo  F(jp  =  Pf)  dpcj  +  pf  j(  F(p  =  pf)  dpcj 


1  rpf  + 1  _  j  TP/  +  1 

+  2  I-  ^  =  pf^ dp^  +  2  -  {Pcj  ~  pf)2  F@  =  pf)  dP*i 

*'P/+l— yro  z  Jpf+1- yro 


N  p2  A2  2  1  i 

k  ~7rA!^5  ~  ^(Pf)}  +  4  Vto [h(Pf)  4-  h)  +  -  Pf[I\  +  h{Pf)) 

~  g  [(exP(~c'2)  ~  (1  ~  Pro)  •  exp(-C2(  1  -  yT0)2)] )  . 


(14a) 


The  right  hand  side  of  Eq.(ll)  becomes, 


rV(p—Pf)  rPQVTQ  ppf 

L,  ..  PdV  =  PfV(p  =  Pf)-[  V{0<p<poyTO)dp+  I  V(pQyTQ<pf)dp\ 

\P  V)  '  JO  J  pnt/Tn  -1 


PqAV  (j/to  +  1) 


Vto  + 1 


(Pf  + 1)  •  h{Pf)  + 


yro  (yro  -  1) 


h-^VTO-h 
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It  should  be  noted  that  the  factors,  N pi  A2s/k,  in  Eq.(14a)  and  pQAV/(yTQ  +  1),  in 
Eq.(14b)  are  identical,  thus  dropping  out  of  the  conservation  of  energy,  Eq.(ll),  as  common 
factor. 

The  p-V  equation  constructed  from  the  mechanistic  model,  Eq.(13b),  contains  five 
unknowns  (A,  po,  AV,  VJy  (or  VL),  yT0),  the  magnitudes  of  which  are  determined  by  min¬ 
imizing  the  differences  between  Eq.(13b)  and  a  specified  data  set  based  on  the  method 
of  least  squares,  under  the  constraint  imposed  by  the  conservation  of  energy,  Eq.(ll)  and 
Eq.(14).  Because  the  p-V  equation  consists  of  two  equations,  one  for  the  high  pressure 
region  and  the  other  for  the  low  pressure  region,  and  also  because  algebraic  equations  re¬ 
sulting  from  the  application  of  the  method  of  least  squares  are  non-linear,  a  computational 

I  - 

program  is  developed  to  find  the  five  unknowns.  The  program  requires  a  set  of  initial  guess 
values  for  the  five  unknowns.  The  parameters,  A,  pQ,  A  V,Vu,  of  the  error  function  p-V 
equation,  Eq.(lb),  are  used  for  initial  values  with  the  initial  value  for  the  fifth  unknown, 
pro,  being  set  to  zero.  The  program  employs  the  Newton-Raphson  iterative  technique 
around  the  value  of  yro  to  minimize  the  errors  between  Eq.(13b)  and  the  data  points 
while  conservation  of  energy  is  satisfied  exactly,  until  the  five  unknown  converge  to  a  set 
of  solutions. 

Discussion  of  Results 

We  begin  with  physical  interpretations  of  parameters  of  p-V  equations  in  terms  of  the 
'mechanistic  model.  The  first  equation  in  Eq.(12)  is, 

AV  =  NV0(yTQ  +  l) 

=  N(V0  +  AsyT). 

Noting  that  yro  Is  a  ratio  of  the  piston  displacement  by  elastic  tissue  distension  to  that 
by  alveolar  recruitment,  AV  of  the  error  function.  p-V  equation  is  the  maximum  possible 


volume  available  for  inflation,  and  is  related  to  the  mechanistic  model  as  a  product  of  the 
total  number  of  elements,  N,  and  the  elemental  volume  available  for  inflation  through  both 
the  ‘pop-open’  mechanism,  V0  (corresponding  to  the  volume  increase  due  to  the  alveolar 
opening),  and  the  piston  displacement,  A3yT  (corresponding  to  the  elastic  wall  distension  of 
TRS).  Vu  is  related  only  to  the  solution  of  the  high  pressure  region  as  Vu  =  V  (p  — >  co). 
On  the  other  hand,  under  the  two-region  p-V  equation  of  the  mechanistic  model,  the 
definition  of  AV  needs  to  be  elaborated.  Since  VL  #  V  (p  — >  —  oo)  in  the  lower  pressure 
solution  of  the  mechanistic  model,  AV  should  be  interpreted  as  the  maximum  possible 
volume  change  when  the  high  pressure  solution  is  extended  into  the  low  pressure  region. 

The  second  equation,  p0  =  ( k/Aa )  y0,  indicates  that  the  pressure  at  the  midpoint  of  the 
p-V  curve  is  an  equivalent  pressure  required  to  displace  the  piston  against  the  spring  force 
over  the  pop-open  displacement  of  y0.  It  may  be  rewritten  as  p0V0  =  kyl,  therefore,  p0V0/ 2 
is  the  spring  energy  required  to  displace  the  piston  by  the  amount,  yQ.  This  observation 
implies  that  the  pressure,  pQ,  is  related  to  both  the  alveolar  recruitment  (through  y0)  and 
the  elastic  tissue  distension  (through  k).  A  higher  magnitude  of  pQ  implies  a  larger  value 

of  the  spring  constant  9wall  elasticity)  and/or  a  greater  amount  of  energy  required  to 
pop-open  the  elements. 

The  non-dimensional  parameter,  A,  is  related  to  the  parameter,  C,  of  the  mechanistic 
model  through  the  third  equation  in  Eq.(12),  C  =  (VSF/4)  A.  As  may  be  seen  from  Eq.(9), 
the  parameter,  C,  appears  as  a  factor  in  the  function,  I3(p).  Since  the  function,  I3(p),  is  a 
monotomcally  increasing  function  of  p,  an  increase  in  volume,  V,  becomes  more  sensitive 
to  a  change  in  pressure  when  the  magnitude  of  A  is  larger.  The  observations  above  may 
be  further  extended  in  terms  of  the  distribution  of  elements  over  the  critical  pop-open 
pressure,  Eq.(13c).  The  number  distribution  of  elements  is  a  normal  distribution  with  its 
mean  at  p(=  p/pQ  -  l)  and  a  standard  deviation,  a,  which  is  proportional  to  1/A.  Since 
the  peak  of  the  distribution  is  located  at  pcj  =  p,  the  rate  of  increase  in  the  number  of 
open  elements  increases  (decreases)  for  p  <  p0  (p  >  p0);  an  observation  consistent  with 
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the  fact  that  p0  is  a  pressure  at  the  inflection  point  of  the  error  function  p-V  equation. 
A.  larger  value  of  A  indicates  a  smaller  standard  deviation,  indicating  a  higher  pealc  in 
number  density  and  a  sharper  distribution. 

The  p-V  equation,  Eq.(13b),  of  the  mechanistic  model  consists  of  the  low  pressure 
solution  in  which  all  open  elements  are  unsaturated,  and  the  high  pressure  solution  where 
some  elements  are  saturated  (fully-distended).  The  equation  has  three  pressure  -dependent 
terms.  A  term  proportional  to  p  in  the  equation  for  the  low  pressure  region  (the  third 
term)  is  due  to  the  elastic  distension  of  the  elements  that  are  open  at  p=0.  Two  other 
pressure-dependent  functions  are  I2(p),  originating  from  volume  changes  due  to  the  piston 
displacement,  and  /3(p),  which  results  from  both  the  pop-open  volume  and  the  piston 
displacement.  The  former  is  symmetric  with  respect  to  p(=  p/p0  -  1)  =  0,  i.e.  I2(p)  = 
h(-p)]  while,  the  latter  is  antisymmetric,  i.e.  /3((p)  =  -  h{-p).  Furthermore,  the  p-V 
equation  in  the  high  pressure  region,  V  {p0-yTo  <  p),  is  independent  of  the  magnitude  of 
2/to  i  while,  the  solution  V),  is  sensitive  to  the  magnitude  of  Uto  in  both  the  low  and  the 
high  pressure  regions. 

Fig-5  shows  six  representative  data  sets  of  patients  with  ARDS  as  well  as  the  cor¬ 
responding  p-V  equation,  Eq.(13b),  derived  from  the  mechanistic  model  for  the  inflation 
(I)  process.  Fig. 6  is  a  plot  of  [the  volume  predicted  by  model-based  p-V  equation  at  a 
specified  pressure]  vs  [the  corresponding  data  volume]  for  all  inflation  data  points  from 
the  twenty  one  data  sets.  Both  figures  show  very  good  agreements  between  the  model  and 
the  clinical  data  with  R~  for  Fig.  6  being  equal  to  0.9993.  (Various  parameters  for  all 
data  sets  are  summarized  in  Table  1.)  The  solid  (dotted)  curves  in  Fig.5  are  the  solution 
of  the  low  (high)  pressure  region  (i.e.  the  first  (second)  equation  in  Eq.(  13b))  with  the 
composite  solution  indicating  that  the  p-V  curve  is  not  antisymmetric  with  respect  to  po- 
However,  since  the  error  minimization  is  applied  between  the  antisymmetric  error  function 
p-V  equation,  Eq.(lb),  and  the  mechanistic  model  equation,  Eq.(13b),  the  two  curves  are 
very  close  to  each  other  in  the  low  pressure  region  of  0  <  p  <  PqVto  •  Ranges  of  various 
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parameters  listed  in  Table  1  are, 

A  =  1.5  -  5.5,  po  =  13  -  31  [cmH20], 

Ay  =  1-4  [L],  yT0  =  0.289  -  0.695. 

Since  the  boundary  pressure  between  the  low-  and  the  high-  pressure  solution,  pQ  ■  yT0,  is 
low  compared  to  the  end-of  inflation  pressure  for  most  data  sets,  the  antisymmetric  high 
pressure  solution  is  applicable  over  a  major  part  of  the  data  sets  analyzed.  It  is  also  noted 
here  that  if  the  condition  of  yT0  —  0  (negligible  elastic  tissue  distension)  is  imposed,  the 
solution  of  the  mechanistic  model,  consisting  solely  of  the  solution  for  the  high  pressure 
region  becomes  identical  to  the  (antisymmetric)  error  function  p-V  equation  although  the 
conservation  of  energy  is  not  satisfied  by  the  condition.  Fig.  7  is  presented  to  show  the 
magnitude  of  the  left  hand  side 'of  conservation  of  energy  divided  by  p0AV  as  the  abscissa, 
and  (-l)-(the  right  hand  side  of  conservation  of  energy  divided  by  PqAV)  as  the  ordinate 
for  all  data  sets  when  i/to  is  set  to  zero  and  the  parameters  (A,  p0,  AV,  V£)  of  the  error 
functidn  p-V  equation  are  used  for  the  evaluation.  The  figure  shows  that  the  left  and  right 
hand  side  of  conservation  of  energy  have  opposite  signs  for  all  data  sets,  indicating  that 
conservation  of  energy  is  not  satisfied  at  yTQ  =  0. 

The  range  of  yT0  obtained  by  the  mechanistic  model  indicates  that  the  fraction  of  total 
volume  available  for  the  pop-open  mechanism  (alveolar  recruitment),  NVq/AV,  which  is 
equal  to  1/(1  +  yro),  ranges  between  0.59  and  0.78. 

Summary 

A  mechanistic  model  of  TRS  elements,  each  consisting  of  a  piston-spring  system,  is 
developed  to  analyze  quasi-static  pressure-volume  curves  for  the  inflation  process.  The 
model  accommodates  both  the  alveolar  recruitment  (in  terms  of  the  critical  pop-open 
pressure)  and  the  elastic  distension  of  wall  tissues  (in  terms  of  the  piston  displacement). 
Model-based  relations  (Eq.(12))  are  established  between  the  parameters  in  the  p  —  V  curve 
represented  by  the  error  function  equation,  Eq.(l),  and  in  the  mechanistic  model.  Under 
the  constraint  imposed  by  conservation  of  energy,  the  parameters  of  the  model-based  p-V 
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equation  is  determined  for  each  clinical  data  set  by  a  computational  minimization  proce¬ 
dure  between  the  equation  and  the  data  points,  results  of  which  show  excellent  agreements 
between  the  two  (Figs.5  and  6).  The  p-V  equation  thus  derived,  Eq.(13b),  consists  of  two 
equations;  one  for  the  low  pressure  region  where  all  open  elements  are  active  (=  unsatu¬ 
rated)  as  the  piston  of  an  element  is  yet  to  reach  its  stroke  limit,  and  the  other  for  the  high 
pressure  region  where  some  open  elements  are  saturated.  The  elemental  distribution  over 
.the  critical  pop-open  pressure,  Eq.(13c),  is  a  normal  distribution  with  its  shape  (the  mean 
and  the  standard  deviation)  affected  substantially  by  the  magnitudes  of  two  parameters 
in  the  mechanistic  model,  A  and  pq. 

The  present  analysis  is  for  the  inflation  process.  The  deflation  process  is  different 
from  the  preceding  inflation  because  of  the  absence  of  the  pop-open  mechanism,  and  also 
becasue  of  a  possibility  of  airway  closure.  However,  a  certain  aspect  of  the  deflation  process 
may  be  predicted  from  the  inflation  analysis;  which  will  be  discussed  in  Prat  II  as  a  validity 
test  of  the  mechanistic  model.  Wide  ranges  covered  by  the  parameters,  A,  po ,  A  V  and 
Uto  ■  of  the  p-V  equation  and  the  mechanistic  model  need  to  be  interpreted  in  terms  of  the 
shape  and  the  range  of  the  p-V  curves  as  well  as  in  terms  of  the  elemental  distribution  and 
its  changes  along  the  corresponding  p-V  curve;  which  will  also  be  discussed  in  Part  II. 
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Figure/Table  Captions 

Table  1.  Summary  of  Inflation  Data  Analyses. 

Fig.l.  A  typical  quasi-static  pulmonary  pressure- volume  curve. 

Fig. 2.  Error-function  p-V  equation  and  inflation  data  points. 

Fig.3.  Error-function  p-V  equation  and  deflation  data  points. 

Fig. 4.  A  schematic  diagram  of  mechanistic  model  of  TRS  element. 

Fig. 5.  Model-based  p-V  equation,  Eq.(13b),  vs  data  points  for  inflation  process, 
solid:  solution  for  low  pressure  region, 
dotted:  solution  for  high  pressure  region. 

Fig. 6.  V  (volume  predicted  by  model-based  p-V  equation)  vs  V  (volume  of  data) 

I  -  • 

for  a  specified  pressure. 

Fig. 7.  AU/po  •  AV  vs  (—1)  ■  f  pdV/po  •  AV  when  yro  =  0 
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Table  1.  Summary  of  Inflation  Data  Analysis. 


Data  No 

A 

A.l. 

2.9652 

2. 

2.9578 

B.l. 

2.7173 

2. 

2.7304 

C.l. 

3.3532 

2. 

3.3664 

D.l. 

1.6273 

2. 

1.9257 

E.l. 

2.7160 

2. 

2.6497 

F.l. 

2.6029 

2. 

2.6423 

G.l. 

1.7288 

2. 

1.9277 

H.l. 

1.8901 

2. 

2.0421 

I.  1. 

3.5379 

2. 

3.5449 

J.  1. 

2.7981 

2. 

2.7364 

K.l.  ' 

2.4296 

2. 

2.4708 

L.  1. 

1.2500 

2. 

1.5318 

Ml. 

3.2725 

2.. 

2.9972 

N.l. 

2.7915 

2. 

2.8046 

0.1. 

1.9487 

2. 

2.1412 

P.l. 

2.4566 

2. 

2.4553 

Q.l. 

2.1027 

2. 

2.2041 

R.l. 

1.1672 

2. 

1.6209 

S.l. 

3.1381  - 

2. 

3.0894 

T.l. 

5.4709 

2. 

5.4708 

U.l. 

3.2818 

2. 

3.2819 

Po  AV 

[cmH20]  [L] 


22.398 

2.3726 

22.411 

.  2.3772 

21.981 

1.5612 

21.999 

1.5559 

25.082 

1.6193 

25.073 

1.6139 

13.324 

3.1567 

13.999 

2.8392 

30.361 

1.6216 

30.817 

1.6847 

23.880 

1.5066 

23.863 

1.4874 

20.156 

3.0989 

20.731 

2.8704 

14.959 

1.7905 

15.405 

1.7129 

25.248 

2.7570 

25.232 

2.7508 

26.208 

3.7326 

26.830 

3.9129 

17.895 

1.3424 

17.951 

1.3304 

11.592 

1.2470 

13.213 

1,1256 

29.865 

3.9075 

30.327 

4.2463 

15.310 

1.6284 

15.297 

1.6219 

18.374 

1.8797 

18.837 

1.7753 

26.982 

1.3277 

27.032 

1.3278 

19.314 

1.3316 

19.583 

1.2941 

13.986 

2.0685 

16.352 

1.7396 

26.802 

3.1306 

26.925 

3.1878 

30.038 

1.7695 

30.037 

1.7694 

24.439 

2.8956 

24.349 

2.8956 

Vto 

-0.0709 

-0.0723 

.  0.347 

-0.0762 

-0.0728 

0.359 

-0.0365 

-0.0342 

0.329 

-0.5727 

-0.3542 

0.379 

-0.0768 

-0.0923 

0.480 

-0.0962 

-0.0873 

0.365 

-0.3645 

-0.2138 

0.447 

-0.2185 

-0.1563 

0.431 

-0.0350 

-0.0333 

0.368 

-0.1583 

-0.1887 

.  0.474 

-0.0997 

-0.0916 

0.440 

-0.2829 

-0.1678 

0.695 

-0.0861 

-0.2179 

0.626 

-0.0705 

-0.0681 

0.358 

-0.2423 

-0.1696 

0.389 

-0.0787 

-0.0773 

0.467 

-0.1279 

-0.1004 

0.389 

-0.5149 

-0.2530 

0.406 

-0.0851 

-0.1030 

.  0.407 

-0.0097 

-0.0097 

0.289 

-0.2075 

-0.2075 

0.327 

1.  Obtained  by  applying  the  method  of  least  squares 
along  with  error  function  p-V  equation. 

2.  Results  from  the  mechanistic  model. 
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& 


Fig. 2.  Error-function  p-V  equation  and  inflation  data  points. 


Fig. 5.  Model  based  p-V  equation,  Eq.  (13b),  vs  data  points  for  inflation  process. 
Solid:  solution  for  low-pressure  region, 

Potted:  solution  for  high-pressure  region. 
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Fig. 6.  V  (volume  predicted  by  model-based  p-V  equation)  vs.  V(Volume  of  data)  for  a  specified  pressure. 
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Fig.7. 
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Chapter  3.  A  Mechanistic  Model  for  Quasi-Static  Pulmonary  Pressure- Volume 
Curves:  Examination  of  Clinical  Data 

Abstract 

A  p-V  equation  is  developed  in  Part  I  based  on  a  mechanistic  model  of  a  total  respiratory 
system.  In  Part  II,  twenty  one  p  —  V  data  sets  of  patients  with  acute  respiratory  distress 
syndrome  are  examined  using  the  mechanistic  model,  relating  the  quasi-static  pulmonary 
p-V  curve  to  the  corresponding  respiratory  conditions  in  terms  of  a  volume  increase  due 
to  alveolar  recruitment  and  due  to  elastic  tissue  distension,  the  elemental  distribution 
ranging  from  the  closed  elements  to  the  saturated  (open  and  fully-distended)  elements 
and  its  changes  with  pressure.  The  compliance  (local  gradient)  of  p-V  curves  is  shown 
to  represent  the  change  in  the  total  volume  of  saturated  elements;  while  the  pressure  at 
the  maximum  compliance  is  identified  as  the  location  where  a  maximum  rate  of  increase 
occurs  both  in  the  volume  increase  due  to  alveolar  recruitment  and  in  the  volume  increase 
due  to  an  increase  in  the  saturated  elements.  Validity  of  the  model  is  provided  by  its 
predictions  of  the  upper  volume  asymptote  and  the  maximum  possible  volume  change  of 
the  corresponding  deflation  process  which  agree  well  with  the  clinical  data. 
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Nomenclature 


A, 

B 

C 


/,  F 

Ii(i  =  1-5) 

k 

N 

Nopen 
Nsat 
NUn  sat 

Nj 

P 

i 

P 


piston  surface  area  on  which  pressure  is  acting 
(k/As)yx  =  p0  -yT0 

j 

=  v^A/4 

distribution  functions  (Eq. (3)) 
functions  defined  in  Eq.(2) 
spring  constant  [N /m] 
total  number  of  TRS  elements 
total  number  of  open  elements 
total  number  of  saturated  elements 
total  number  of  unsaturated  elements 
number  of  elements  at  energy  level  j 
pressure  (interpleural  pressure  difference) 
non-dimensional  pressure,  p/p0  -  1 


Pcj 


critical  pressure  at  which  an  element,  j,  ‘pops  open’. 


Pcj 

Pf 

Vo 


PID 

V 


'pop— open 

Vsat 

Vopen—sat  (p,  =  0) 


Vn 


V, 


L(U) 


vS 


Pcj/Po 

pressure  at  the  end  of  inflation 

pressure  at  the  inflection  point  in  model  equation 

pressure  at  the  intersect  of  inflation  and  deflation  processes 
volume 

total  ‘pop-open  volume’ 

total  volume  of  saturated  elements 

total  volume  of  elements  open  at  p  =  0  when  they  are  all  saturated, 
volume  change  from  the  state  of  p  =  0 
lower  (upper)  bound  of  volume 

an  (imaginary)  upper  bound  of  volume  for  the  deflation  process 
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V 

non-dimensional  volume,  ( V  -  ((Vu  +  Vl)/2)/(AV/2) 

Vj 

volume  of  an  element  j 

AV 

Vu-VL  =  NVo{yTo  +  1) 

% 

‘pop-open’  volume  (=  Agy0 ) 

Vj 

piston  displacement  of  an  element  j 

>-S,  \ 

yo 

=  Vo /As 

Vt 

piston  stroke  limit 

Vto 

Vt/vq 

Greek  symbols: 

a 

constant  of  proportionality 

A 

apo AV (non-dimensional  parameter)  (Eq.(lb)) 

a 

(8/ 7r) 1/2/ A,  Standard  deviation  (Eq. (3b)) 

&D 

cr -po 

Superscript: 

- 

— 

related  to  a  single  TRS  element 

d . 

deflation  process 

Acronym: 

ARDS 

acute  respiratory  distress  syndrome 

TRS 

total  respiratory  system 

(Eq.(lb)) 
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Introduction 


In  Part  I,  the  error  function  p-V  model  equation  is  shown  to  agree  well  with  clinical 
p-V  data.  The  equation  and  the  corresponding  non-dimensional  form  are, 

7=erf(Cp).  (la, 6) 

where  AV  =  Vy  -  VL,  Vy  —  upper  volume  asymptote,  VL  =  lower  volume  asymptote, 

Po  —  a  pressure  at  the  inflection  point  of  the  curve,  C  =  y/nA/i, 

V  =  [V  -  (Vu  +  Vt)/2]/(AV/2),  p  =  (p/p0)  -  1. 

The  mechanistic  model  of  a  TR.S  element  developed  in  Part  I  is  a  piston-spring- 
cylinder  system  with  the  alveolar  recruitment  and  the  elastic  tissue  distension  represented 
respectively  by  the  critical  pop-open  pressure,  pCj,  and  by  the  displacement  of  piston 
against  the  spring  force,  yj  =  A2ip—pCj)/k  =  piston  surface  area,  k  =  spring  constant). 

Based  on  the  error  function  p-V  equation  and  the  mechanistic  model,  and  allowing 
for  a  distribution  of  elements  over 'pcj,  the  following  model-based  p-V  equations  as  well  as 
the  corresponding  distribution  function  of  TRS  elements  are  derived, 


V  ( 0  <  p  <  Po  •  yro )  = 


Vu  +  VL 


+  +  + 1  -  M  ^ 

t/  /  ~  '  >  Vu  +  VL  AV 

v  (po-Vto  <p)  = - - - +  —  hip), 


,  (2a) 


i — txp(- 

7T  <?d  \  2  L  crD  J  J 


N-dp 


where 


h  =  erf(C),  I2{p)  =  expi~C2 p2)  -  expi~C2),  /3(p)  =  erfiCp), 
h  =  erfiCil  -  yTQ)),  J5  =  exp(-C2  (yT o  -  l)2)  -  exp(-  C2), 

°D  =  (8/tt)2  p0/A,  cr  =  (8/tt)5/A  pcj  =  pCJ/p0. 
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Following  the  development  of  the  mechanistic  TRS  model  and  the  discussion  on  relations 
between  the  error  function  p-V  equation  and  the  the  model-based  p-V  equation  as  well  as 
between  parameters  of  the  two  equations  in  Part  I,  our  discussions  here  concentrate  on 


various  results  that  may  be  derived  from  applications  of  the  model  to  the  clinical  p-V  data 
sets. 

Relationship  between  Inflation  and  Deflation  Processes 

Although  this  report  is  focused  on  the  mechanistic  model  for  the  inflation  process, 
there  exist  certain  relations  between  the  inflation  and  the  deflation  process  that  may  be 
evaluated  from  the  present  inflation  analyses.  We  consider  a  general  case  in  which  a 
quasi-static  inflation  process  proceeds  to  a  pressure,  (=  end-of-inflation  pressure  = 
initial  pressure  of  the  corresponding  deflation  process),  followed  by  a  quasi-static  deflation 
process.  In  terms  of  the  mechanistic  model,  TRS  elements  at  p  =  Pid  with  its  critical 
pop-open  pressure  less  than  zero  ( p_Cj  <  0)  are  still  closed  and  have  not  contributed  to 
the  volume  change  during  the  inflation  process  from  p  =  0.  to  p  =  Pid\  hence,  we  may 


postulate  that  only  those  elements  that  axe  open  at  p  =  pjp  participate  in  the  deflation 
process  to  follow.  Therefore,  Vfi  (=an  (imaginary)  upper  bound  of  volume  for  the  deflation 
p-V  curve)  may  be  viewed  as  the  volume  which  would  be  attained  if  the  elements  that  are 
open  at  the  end  of  the  inflation  process,  p  =  piDj  were  all  fully  saturated;  i.e. 


vS  =  v(p  =  p,D)  + 


roo 

/  Wo(l  +  vto)  ■  F(p  =  pID)  dpCj  -  vp(p  =  pID) 
Jo 


The  first  term  on  the  right  hand  side  is  the  inflation  volume  at  p  —  Pid*  The  second  integral 
term  is  the  volume  summed  over  all  open  elements  at  p  —  pid  when  they  are  saturated, 
and  the  last  term  is  the  actual  volume  increase  in  the  inflation  process  from  p  =  0  to 
V  —  Pi d  with  the  two  terms  in  the  square  bracket  together  representing  a  volume  increase 
above  V(p  =  pID)  if  all  open  elements  at  pID  were  saturated.  Under  the  assumption 
that  the  magnitude  of  pjo  is  greater  than  B  (=  p0  yT0)  which  is  valid  for  all  data  sets 
analyzed,  Eqs.(2,3)  along  with  Eq.(9)  of  Part  I  for  Vp  are  used  to  evaluate  the  right  hand 
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side,  yielding 


Yu  =  Vl  + 


+ 


2AV 


7rA(l  +  Utq) 
AV 


+  ~T~  •  t1  +  h(Pw)l 


2(1  +  yT0)  t1  ~  (2  ~  Vto)  •  h  +  (1  -  m)  •  It] 


(4) 


The  clinical  data  sets  made  available  to  us  contain  both  the  inflation  and  the  deflation  p-V 
curves  for  each  patient  with  ARDS;  however,  the  p-V  curves  are  obtained  separately  for  the 
inflation  and  the  deflation  process.  (See  [1]  for  the  procedure  of  data  acquisition.)  Fig.l 
shows  inflation  (unfilled)  and  deflation  (filled)  data  points,  as  well  as  the  corresponding 
inflation  (I)  and  deflation  (D)  curves  for  a  typical  data  set  we  examined.  The  inflation 
curve  in  Fig.l  is  Eq.(2)  of  the  mechanistic  model;  while,  the  deflation  curve  is  obtained 
by  straight  applications  of  the  method  of  least  squares  between  data  points  and  the  error 
function  p-V  equation,  Eq.(l).  As  may  be  observed  from  Fig.l,  the  end-of-inflation  point  is 
quite  different  from  the  initial  deflation  point  for  most  data  sets.  To  accommodate  the  data 
into  the  analysis  based  on  Eq.(4),  the  initial  deflation  data  point  is  translated  horizontally 
until  it  meets  the  inflation  p-V  curve,  the  pressure  value  of  which  is  then  defined  as  pjd 
in  Eq. (4) ,  as  indicated  in  Fig.l,  implying  that  the  deflation  curve  preceded  by  an  inflation 


curve  is  assumed  to  be  the  same  as  the  deflation  curve  of  data  sets  horizontally  translated 
until  the  beginning-of-deflation  data  point  is  on  the  inflation  curve. 

Fig. 2  presents  V$  of  Eq.(4),  predicted  from  the  mechanistic  model  of  the  inflation 
process,  plotted  against  Vtf  of  the  error-function  p-V  equation,  Eq.(l).  (For  a  complete 
list  of  numerical  results  relevant  to  the  analysis,  see  Table  1.)  A  maximum  and  a  minimum 
of  a  difference,  Vj}  (Eq.(4))  -  V*  (Eq.(  1)),  are  0.1113  [L]  and  -0.0352  [L]  respectively  with 
an  average  of  the  difference  =  0.0460  [Lj.  Agreements  are  very  good  in  view  of  the  fact 
that  Eq.(4)  predicts  the  upper  volume  asymptote  of  the  deflation  process  in  terms  of  the 
conditions  predicted  by  the  mechanistic  model  of  the  corresponding  inflation  process;  thus 
indirectly  supporting  a  certain  degree  of  validity  of  the  mechanistic  model.  Also,  the  fact 
that  the  magnitude  of  p/£>  is  determined  from  the  horizontal  translation  of  the  deflation 
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curve  indicates  that  the  deflation  process  may  be  relatively  insensitive  to  the  inflation 
history  prior  to  p  =  pID.  On  the  other  hand,  Eq.(4)  for  V$  contains  pID  as  variable  for 
specified  inflation  conditions,  indicating  that  the  magnitude  of  V$  (i.e.  the  shape  of  the 
deflation  curve)  changes  as  the  end-of-inflation  pressure  is  varied.  A  similar  statement 
has  been  made  previously  by  Jonson  [2].  A  closer  examination  of  Fig.2  indicates  that 
the  mechanistic  model  slightly  underpredicts  V$  compared  to  that  of  the  error-function 
equation  for  most  of  data  sets.  This  could  indicate  either  a  quantitative  limitation  of  the 
mechanistic  model  or  the  effects  of  the  inflation  process  preceding  the  deflation.  If  yT0  is 
set  to  zero,  Eq.(4)  is  reduced  to 


VS  @to  =  0)  =  Vu  -  AV  (/,  -  h(p1D))/ 2.  (4a) 


Fig.3  plots  VL1  of  Eq.(4a)  vs  of  the  error  function  p-V  equation  for  the  deflation  process. 
Agreements  are  fairly  good  between  the  two.  Results  presented  in  Fig.3  reflects  that  the 


p-V  curve  is  relatively  insensitive  to  the  magnitude  of  Uto  i  a  reason  why  the  antisymmetric 
error  function  p-V  equation  (for  which  yT0  is  zero)  fits  well  with  p-V  curves. 

The  error  function  p-V  equation  fits  well  not  only  with  the  inflation  but  also  with  the 
corresponding  deflation  processes  as  shown  in  Part  I.  Therefore,  regardless  of  the  actual  de¬ 
flation  process,  we  may  define  AVd  ( A  V  of  the  deflation  process)  as  the  maximum  possible 
volume  change  of  a  specified  TRS  during  the  deflation  process;  which,  in  our  mechanistic 
model,  yields  the  following  equation  for  AVd  in  terms  of  the  inflation  parameters: 


AVd  =  V0  (yr o  +  1)  f  N  ■  F(p  =  pID)  dpcj 
Jo 


AV 


(1  +  h(pi£))). 


(5) 


Fig.4  is  AV  of  deflation,  AVd,  predicted  by  Eq.(5)  plotted  against  the  corresponding  AVd 
of  error-function  p-V  equation  determined  by  the  method  of  least  squares.  (See  Table  1 
for  numerical  values.)  Agreements  between  the  two  are  reasonably  good  for  a  majority 
of  data  sets,  except  for  six  data  sets  shown  in  filled  circles  accompanied  by  alphabetical 
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data  numbers.  Shown  in  Fig.  5  for  the  six  data  sets  are  their  deflation  data  points  and 
two  p-V  curves;  one  (dotted)  for  the  error  function  p-V  equation  with  the  method  of  least 
squares  applied  to  determine  (Ad,  p$,  &Vd,  V$),  and  the  other  (solid)  for  the  error  function 
p-V  equation  with  (AFd,  Vjf)  determined  from  Eqs.(4,5)  and  (Ad,  p$)  determined  by  the 
method  of  least  squares.  Two  curves  are  different  in  their  approaches  to  different  high  and 
low  asymptotes.  Although  agreements  of  the  solid  curves  (with  two  adjusting  parameters) 
with  the  data  points  are  not  as  good  as  that  of  the  dotted  curves  (with  four  adjusting 
parameters),  the  errors  are  small  for  the  solid  curves  in  view  of  the  fact  that  the  magnitude 
of  AFd  is  quite. different  between  the  two  curves,  indicating  that  a  better  understanding  of 
relations  between  the  p-V  equation  and  the  corresponding  intra-respiratory  changes  helps 
interpret  various  characterisitcs'of  p-V  curves  accurately. 

Interpretation  of  Inflation  p-V  Curves  based  on  Mechanistic  Model 

Fig. 6  shows  ranges  covered  by  all  data  sets  analyzed  in  terms  of  po  (the  inflection 
pressure  of  the  high-pressure  solution),  AF  (maximum  volume  available  for  inflation  in 
the  high-pressure  solution)  and  1/ (1  -t-yro)  (the  fraction  of  total  elemental  volume  available 
for  the  pop-open  mechanism  (alveolar  recruitment),  NV0/AV),  all  plotted  against  the  non- 
dimensional  parameter,  A.  The  data  sets  with  their  alphabetical  data  numbers  indicated 
in  the  figure  are  those  to  be  analyzed  in  detail  in  comparative  analyses  to  follow.  (Various 
parameters  of  the  six  data  sets  are  reproduced  as  Table  2.  Parameters  of  all  data  sets 
are  listed  in  Table  1  of  Pant  I.)  The  range  of  A  is  between  1.5  and  3.6  except  for  Data 
T  (A  =  5.47).  The  six  data  sets  (B,  E,  M,  N,  R,  T)  are  different  from  each  other  in  the 
following  ways: 

1.  Data  set  B,  E  and  N  have  roughly  the  same  magnitude  in  both  A  (  2.65  -  2.80)  and 

AF  (=~  1.55  -  1.68)  with  significantly  different  values  for  p0. 

2.  Data  set  E  and  T  have  substantially  different  values  of  A  with  po  and  AF  being 
approximately  the  same  in  magnitude. 
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3.  N  and  R  also  show  similar  characteristics  of  being  different  in  A  and  common  in  p0  and 
AV ;  however,  the  pair  covers  lower  range  in  both  A  and  po  compared  to  E  and  T. 

4.  E  and  M  are  different  in  terms  of  the  magnitude  of  AV  with  p0 e  ~  Pom,  A e  Am- 

5.  T  and  R  represent  the  data  sets  with  very  high  and  low  values  of  A,  respectively. 

Figs. 7  and  8  are  p-V  and  the  corresponding  non-dimensional  p-V  curves  over  their  ranges 
of  measurement  for  Data  Set  B,  E  and  N  in  (a),  Data  Set  E  and  T  in  (b),  Data  Set  N 
and  R  in  (c)  and  Data  Set  E  and  M  in  (d). 

Referring  to  Fig.7  and  noting  that  the  dotted  (broken)  vertical  lines  are  the  location 
of  p  =  po  (po  •  Pro),  it  may  be  observed  that  the  range  for  the  low  pressure  solution, 
0  <  p  <  po  •  pro,  in  which  all  elements  are  active. (unsaturated),  is  very  narrow  compared 
to  the  range  for  the  high  pressure  solution  with  an  exception  of  Data  set  M,  for  which  the 
measurement  does  not  reach  the  inflection  pressure,  po,  with  p0  •  pro  >  15  [cmH20\. 

The  corresponding  non-dimensional  (p  -  V )  curves,  based  on  the  definitions  of  V  and  p  in 
Eq.(l),  represent  Eqs.(2a,b)  in  the  following  normalized  form: 

V  (—1  <p  <  yT0  —  1 )  = 

+  (/a®  +  ^ 

^(Pro-l<p)  =  h(P),  (6) 

The  normalization  of  volume  transforms  the  two  volume  asymptotes,  V\j  and  Vl  into 
V  =  +1  and  V  =  -1  respectively;  while,  the  pressure,  p  =  p0,  at  the  inflection  point 
is  transformed  into  p  =  0.  With  both  the  location  of  po  and  the  volume  asymptotes 
made  common  to  all  p-V  curves,  the  resulting  non-dimensional  representations  in  Fig. 8 
are  characterized  by  a  single  non-dimensional  parameter,  A.  The  parameter,  A,  is  twice 
the  maximum  local  compliance  at  p  =  po  (dV / dp(p  =  0)  =  A/2).  Since  the  compliance 
is  maximum  at  the  origin  of  p  —  V  diagram,  the  first  quadrant  ( V ,  p  >  0)  in  Fig. 8  is  a 
region  of  decreasing  local  compliance  with  pressure;  while,  the  third  quadrant  (V,  p  <  0) 
is  a  region  of  increasing  local  compliance  with  pressure.  The  origin  (p  =  0,  V  =  0)  of 


^QM-1-  yT0)-™f±I4  +  2-A 


7rA 
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dimensional  p-V  curves  is  transformed  into  (p  =  —1,  V(V  =  0))  on  a  p-V  curve;  hence, 
the  physiological  lower  limit  of  p  is  —1. 

Fig. 8  (a)  compares  the  three  data  sets,  B,  E  and  N,  among  which  the  magnitude  of  po  is 
substantially  different  with  A  and  A  V  being  approximately  the  same  in  magnitude.  Three 
curves  are  very  close  to  each  other  because  values  of  A  are  similar,  and  the  difference 
between  the  three  appears  as  the  extent  to  which  the  p-V  curves  are  measured  in  the 
region  of  decreasing  compliance  with  pressure.  In  Fig.8  (b)  and  (c)  differences  between 
the  two  data  sets  occur  in  the  magnitude  of  A  (A r  >  Ae,  A/v  >  A R),  resulting  in  the  T- 
and  N-  curves  above  the  E-  and  R-  curves  respectively  in  the  first  quadrant.  Since  the 
magnitude  of  po  for  the  data  set  T  and  E  are  very  high  compared  to  those  for  N  and  R, 
the  region  of  decreasing  compliance  covered  by  the  T-  and  E-  curves  are  narrower  than  N 
and  R.  Because  both  po  and  AV  are  similar  in  magnitudes  between  the  two  data  sets  in 
Fig. 8(b)  and  (c)  the  shape  of  p  —  V  curves  is  very  similar  to  the  corresponding  p-V  curves. 
In  Fig.8  (d)  two  data  sets  with  a  high  value  of  p0  (~  30  [cmH20])  are  shown.  For  the  data 
set  M  the  high  value  of  p0  combined  with  a  high  value  of  AV  limit  the  measured  range 
of  the  p-V  curve  in  the  region  of  increasing  compliance  only,  compared  to  the  data  set  E 
with  a  smaller  magnitude  for  AV. 

Although  the  p  —  V  diagram  helps  distinguish  differences  among  p-V  curves  and 
effectively  bring  out  various  characteristics  of  each  p-V  curve,  it  is  the  information  from 
the  elemental  distribution  that  relates  various  parameters  of  p-V  curves  to  TRS  conditions 
quantitatively.  On  the  normalized  p  —  V  plane  of  Fig.8  the  local  compliance  at  p  =  0 
(dV /dp  (p  =  0))  increases  with  A;  while,  as  the  standard  deviation,  a,  is  proportional  to  l' /A 
in  the  normalized  number  distribution,  Eq.(3b),  the  distribution  becomes  sharper  and  has 
a  higher  peak  as  A  is  increased.  Fig.9  is  a  plot  of  the  number  distribution  (not  normalized) 
vs  the  critical  pop-open  pressure,  pcj  in  [cmi^O],  for  the  six  data  sets  analyzed  in  Figs. 7 
and  8.  The  number  distribution,  dNj/N  ■  dpc  j,  is  a  fraction  in  the  number  of  elements,  the 
critical  pop-open  pressure  of  which  ranges  between  pCj  and  pCj  +  dpCj  in  [1  /cmH20\.  The 
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corresponding  equation  is  Eq.(3a),  which  indicates  that  the  maximum  number  of  elements 
are  present  at  pcj  =  p  -  p0  with  its  magnitude  equal  to  1/V2n  ■  aD  (=  A/4  •  p0).  As  aD 
decreases  the  distribution  becomes  sharper  (cx  l/aD)  and  its  peak  value  (a  aD )  larger. 
(See  Table  2  for  the  magnitude  of  for  each  data  set.)  Two  distributions  are  shown  in 
Fig. 9  for  each  data  set,  one  at  p  —  pq  (with  its  peak  at  p  =  0)  and  the  other  at  p  —  Pf  (end- 
of-inflation  pressure)  (with  its  peak  at  p  =  pf  -  pQ).  The  vertical  line,  p  =  pQ  (pf)  for  the 
distribution  at  p  =  p0  (pf)  indicates  the  pressure  above  which  the  distribution  is  truncated. 
The  dotted  parts  of  the  curves,  as  discussed  in  Part  I,  correspond  to  the  elements  that 
are  open  at  p  =  0.  The  normal  distribution  truncated  at  pcj  =  0  and  pCj  =  p  (i.e.  the 
solid  part  of  the  curves  in  Fig. 9)  shifts  to  the  right  with  an  increase  in  p  as  more  elements 
become  open.  When  p  <  po,  the  peak  of  the  distribution  lies  in  the  negative  range  of  pCj- 
It  should  be  noted  that  in  Data  set  M  (Fig.9  (c))  the  distribution  at  p  =  pf  lies  below  that 
at  p  =  pq  because  the  measured  range  never  reached  p  =  p0. 

An  integral  of  the  distribution  function  over  the  critical  pop-open  pressure  in  Fig.9  should 
yield  various  fractions  in  number  of  elements  (depending  on  the  upper  and  lower  limits  of 
the  integral)  at  each  quasi-static  state.  Also  an  integral  of  a  product  of  the  distribution 
function  and  the  elemental  volume  over  the  critical  pop-open  pressure  should  provide  us 
with  such  quantities  as  the  volume  change  due  to  alveolar  recruitment,  due  to  elastic  wall 
distension  and  due  to  an  increase  in  the  saturated  elements.  The  following  equations  may 
be  obtained  for  changes  in  the  number  fractions: 

Fraction  of  the  number  of  open  elements  at  p  =  p:  Nopen(p  =  p)/N 


(rp/po  \ 

=  JQ  F(P)  dPcj  J  =  [Ji  +  h(p)]/ 2,  /or  p  >  0, 


vfcj  (~  r'*  ^  ^ 

N 


(7a) 


Fraction  of  the  number  of  saturated  elements  at  p  =  p:  Nsat.(p  =  p)/N 


^satXp  ~  p) 
N 


_  rP/Po  yro  \  =  f  [/4  +  h(p)}/ 2,  if  pQ  •  yT0  <  p, 

Jo  \  0,  if  0  <  p  <  p0  •  Vto-  [  } 
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Fraction  of  the  number  of  unsaturated  elements  at  p  =  p:  Nunsat_{p  =  p)/N 
^unaatXp  =  P)  ( -f^open  ^ sat.  A  _  (  [l"i  —  I4]/ 2,  if  Pq  '  PTO  5:  Vi 

N  y  N  )  \  [A  +  J3(p)j/2,  if  0  <p  <p0  -yTQ. 

Fraction  of  the  number  of  open  elements  at  p  =  0:  Nopen(p  =  0 )/N 

jT  F® = -«  -  [1  -  h\n 

It  should  be  noted  that  (1.)  the  equation  for  Nopen(p)/N  does  not  include  the  elements 
that  are  open  at  p  —  0.  (2.)  the  equation  for  Nsat.(p)/N  does  not  account  for  the  elements 
that  are  open  at  p  =  0  and  saturated  subsequently.  (3.)  Nopen(p  =  0)/N  is  a  function  of 
a  single  parameter,  A.  (4.)  the  number  fraction  of  unsaturated  elements,  Nunsat,{p)/N,  is 
independent  of  pressure  for  po  •  Pro  <  p,  indicating  that  as  more  elements  are  recruited  in 
the  region  of  the  high  pressure  solution  the  same  number  of  elements  are  saturated. 
Similarly,  the  following  equations  are  for  volume  changes  as  pressure  is  varied: 

Total  pop-open,  volume:  Vpop-open(p  =  p) 

Vpop-open(p  =  p)  NV\ o  •  ^  =  — —  - •  [ii  4-  I3 (p)]  for  p  >  0.  (8a) 

Total  volume  of  saturated  elements:  Vaat.  (p  =  p) 

V,«.(s>  =  p)  (  E  WVo(l  +  fro)  • 

AV  •  [/4  +  /3(p)]/2,  if  popro  <P,  fi 

0,  if  0  <  p  <  po  •  pro- 

Total  volume  of  elements  open  at  p  =  0  when  they  are  all  saturated:  Vope.n-Sat.  (p  =  0) 

Kp«n-,«.(P  =  0)  (  =  NV0yn  ■  =  °) 

V*  V'TCl 

=  Yrr;r  (1  -  h)  for  p0pro  <  p.  (8c) 

4  1  +  pro 

A  number  fraction  of  open  elements,  Nopentotal  (p  =  p)/N ,  may  be  derived  from  Eqs.(7a) 
and  (7d)  as, 


(7c) 

(7  i) 


N, 


open  total  {jP  p)  f Ngpenip  —  0)  A open {P  —  P)\  1  H"  I3 (p) 


N 


N 


N 


(9) 
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The  function,  I3(p)  defined  in  Eq.(2),  is  an  error  function,  the  magnitude  of  which  depends 
on  A  and  p0.  Therefore,  the  number  fraction  of  all  open  elements,  Nopen  total  ( p  =  p)/N, 
has  the  point  of  antisymmetry  at  p  =  p0  with  two  asymptotes  of  0  and  1.  Fig.10  shows 
Nopen  total  (p  =  p)/N  vs  p  of  the  six  representative  data  sets  in  solid  lines.  Due  to  the 
antisymmetry  with  respect  to  p  =  p0,  the  fraction,  Nopentotal  (p  =  p)/N,  is  exactly  equal 
to  0.5  when  p  is  equal  to  po  as  the  peak  of  the  distribution  is  located  at  p  =  0.  (See 
Fig.9.)  Since  the  location  of  p  =  poPro  is  the  boundary  between  the  low-  and  high- 
pressure  solutions  of  the  p-V  equation,  all  open  elements  are  still  active  and  unsaturated 
for  p  <  PoVto-  The  fraction  of  open  elements  at  p  —  poPro  is  less  than  0.2  for  four 
data  sets  other  than  Data  sets  R  and  M  for  which  the  fraction  is  ~  0.25.  As  pressure 
increases  beyond  PoPto>  some  of  the  open  elements  begin  to  be  saturated,  the  fraction 
of  which,  NaatXp  =  p)/N  of  Eq.(7b)  is  plotted  in  Fig.10  in  broken  lines.  The  number 
fraction  saturated  depends  on  yT o  in  addition  to  A  and  p0.  The  difference  between  the 
two  fractions  plotted  in  Fig.10  varies  with  the  magnitudes  of  the  three  parameters  with 
{Nopentotail N  —  Nsat./N)  at  a  specified  pressure  ranging  from  ~  0.28  for  Data  R  to  less 
than  0.01  for  Data  T.  It  should  also  be  noted  that  the  two  curves  axe  parallel,  indicating 
that  the  rate  of  increase  in  the  number  of  saturated  elements  is  equal  to  the  corresponding 
rate  of  opening  elements,  once  pressure  exceeds  poyro- 

The  number  fraction  of  open  elements  at  p  =  0,  Nopen{p  =  0)/iV,  of  Eq.(7d)  rep¬ 
resents  the  percentage  of  elements  that  only  experience  elastic  displacement,  which  may 
be  interpreted  as  the  elemental  fraction  representing  a  non-alveolar  part  of  TRS  such  as 
airway  tissues  as  well  as  a  dysfunctional  alveolar  part  which  does  not  respond  to  the  re¬ 
cruitment.  Fig.  11  depicts  the  fraction  as  a  function  of  A,  the  only  parameter  affecting  the 
fraction.  As  Eq.(7d)  shows,  the  fraction,  Nopen(p  —  0 )/N,  has  two  asymptotes  of  0  (as 
A  — ►  oo)  and  0.5  (as  A  — *  0).  Fig.ll  indicates  that  the  number  fraction  of  open  elements 
at  p=0  is  very  sensitive  to  the  magnitude  of  A  as  its  value  drops  from  0.5  to  ~  0.1  when  A 
is  changed  from  0  to  2.  It  should  also  be  mentioned  here  that  the  total  volume  of  elements 
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open  at  p  =  0  when  they  are-all  saturated,  Vopen_sat.(p  =  0)  of  Eq.(8c),  depends  on  three 
parameters,  A,  AV  and  yro,  of  the  mechanistic  model.  The  last  column  of  Table  2  hsts 
Vopen—sat.  (p  =  0)  for  the  six  data  sets.  Although  their  magnitudes  are  small  (0.08  [L]  or 
less),  the  value  varies  substantially  among  the  data  sets. 

As  mentioned  in  Part  I,  the  optimization  of  the  ventilator  strategy  is  required  for 
patients  with  acute  lung  injury  in  intensive  care  units  in  terms  of  pressure-  and  volume- 
ranges  [3  -  5] ,  taking  into  account  such  considerations  as  a  change  in  the  recruited  volume 
with  pressure,  a  rate  of  recruitment  and  over  distension  of  the  respiratory  tissues.  Results 
of  the  mechanistic  model  analyses  relevant  to  the  respiratory  ventilation  are  presented  in 
Figs. 12  and  13.  The  volume,  Vpop-open  of  Eq. (8a),  represents  a  volume  change  of  TRS 
due  to  the  pop-open  mechanism-  (alveolar  recruitment)  only  (that  is,  excluding  the  volume 
change  due  to  the  displacement  of  piston  (elastic  tissue  distention));  on  the  other  hand, 
the  volume,  V3at_  of  Eq.(8b),  is  the  total  volume  of  saturated  elements  of  TRS.  They 
both  increase  as  pressure  increases  along  an  inflation  path.  Sketched  in  Fig.  12  are  p  -  V 
curves  (solid,  Eq.(2)),  p  Vp@p— open  curves  (dotted,  Eq.(8a))  and  p  V3Qt  curves  (broken, 
Eq.(8b))  for  the  six  data  sets.  The  vertical  broken  line  in  the  figure  is  the  location  of 
P  —  Po-  The  intersect  between  the  p  —  Vsat  curve  and  the  x  axis  in  the  figure  is  the  pressure 
at  the  boundary  (=  PoJ/to)  between  the  low  pressure  p-V  solution  (in  which  all  open 
elements  are  active  (unsaturated))  and  the  high  pressure  solution  (in  which  a  part  of  open 
elements  are  saturated).  Since  both  the  p  -  Vpop-open  and  the  p  -  V3at  relations  as  well  as 
p-V  (p  >  po  yro)  equation  are  represented  by  a  common  function,  I3  (p)  (=  erf(Cp), 
C  =  VtFA/ 4,  p  =  p/po  —  1),  the  pressure  at  the  inflection  point,  poi  is  not  only  the  pressure 
at  which  the  local  compliance,  dV/dp,  is  maximum,  but  also  the  pressure  location  for  a 
maximum  rate  of  increase  in  Vpop-open  as  well  as  in  Vsat  with  their  gradients  given  by  the 
following  equations; 


dVsatip)  /  dV (jp)  \ 

dp  V  dp  J 


AVA 

4p0 


•  exp(-  ( Cp )2) 


for  poyro  <  P , 
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dVpop— open 


AV  A 


exp(-  ( Cp )2)  for  p  >  0 . 


dp 


1  +  2/to  4po 


(10) 


The  identity  between  dVaat/dp  and  dV/dp  implies  that  the  shape  of  the  p-V  curves  (in  the 
range  of  the  high  pressure  solution)  closely  represents  the  change  in  the  saturated  volume 
rather  than  the  pop-open  volume.  Since  the  gradient  of  Vpop-open  is  smaller  than  that  of 
Vsat  by  a  factor  of  1/(1  +  yro)  (<  1),  the  magnitude  of  Vaat  eventually  becomes  greater 
than  that  of  Vpop-open  as  pressure  increases.  Also,  1/(1  +  yro),  a  fraction  of  the  pop-open 
volume  (=  pop-open  volume/total  volume  of  a  single  TRS  element),  may  be  interpreted 
as  a  gradient  ratio  of  Vpop-open  to  Vsat.  A  smaller  value  of  yro  (he.  smaller  piston  stroke 
limit)  means  that  the  element,  once  it  pops  open,  reaches  the  saturated  state  earlier; 
hence,  for  Data  T  of  Fig.l2(f)  (|/ro  =  0.289,  poyro  =  8.68 [cmH^O])  Vsat  becomes  greater 
than  Vpop~opeTL  at  a  pressure  close  to  po  yro>  while,  for  Data  M  of  Fig.  12(c)  (yro  =  0.626, 
Poyro  =  18.98  [cmH20])  Vsat  does  not  overtake  Vpop-open  within  the  measured  pressure 
range,  i  The  gradients,  dVpop-open/ dp  between  p  =  0  and  p  —  Pf  (=final  pressure)  (solid), 
and  dV sat /dp  between  p  =  po  yro  andp  =  p/  (broken)  are  plotted  in  Fig.  13  for  the  six  data 
sets.  The  gradients  are  symmetric  with  respect  to  po-  The  data  sets  with  high  gradients 
(Data  set  M,  N,  T)  show  high  sensitivity  of  the  gradients  to  pressure  change  near  po.  Other 
data  sets  with  low  gradients,  particularly  Data  set  E  and  R,  indicate  that  the  gradients 
(i.e.  the  local  compliance)  do  not  change  too  much  over  a  substantial  range  in  pressure 
around  po .  , 

Summary 

The  mechanistic  model  of  TRS  developed  in  Part  I  is  applied  to  examine  p-V  curves 
(in  a  form  of  the  error  function  p-V  equation)  of  patients  with  ARDS  with  the  following 
results: 

1.  Parameters  of  the  deflation  process,  V/}  in  Eq.(4)  and  AVd  in  Eq.(5),  predicted  by 
the  mechanistic  model  of  the  corresponding  inflation  process  agree  well  with  those  of 
the  error  function  p-V  equation  for  the  deflation  process  (Figs. 2,  4),  indicating  that  the 
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mechanistic  model  has  a  certain  validity  to  be  used  for  improving  our  quantitative  un¬ 
derstanding  of  various  intra-respiratory  conditions,  that  the  shape  and  characteristics,  of 
deflation  curves  depend  not  so  much  on  the  inflation  history  but  on  the  end-of-inflation 
(the  onset-of-defiation)  pressure,  and  that  relations  between  parameters  in  p-V  equation 
and  TRS  conditions  are  needed  to  strengthen  our  applications  of  p-V  curves  in  clinical 
settings  (Fig.5). 

2.  The  non-dimensional  p-V  curve,  p  —  V  curve,  is  effective  in  distinguishing  differences  in 
magnitudes  of  model  parameters  among  different  p-V  curves  (Figs.7,  8). 

3.  In  the  mechanistic  model,  the  distribution  function,  Eq.(3),  and  its  change  with  pressure 
are  the  basis  for  evaluating  alveolar  recruitment  and  the  elastic  tissue  distension.  The  shape 
of  the  distribution  function  (the  peak  value  and  the  standard  deviation)  is  determined  by 
the  magnitude  of  the  non-dimensional  parameter,  A;  while,  the  magnitude  of  the  pressure 
at  the  maximum  compliance,  po,  and  its  location  relative  to  the  range  of  the  p-V  curve 
are  the  important  factors  affecting  changes  of  the  distribution  with  pressure  (Eqs.(7)(9), 
Figs.9,10,11). 

4.  In  addition  to  A  and  p0,  other  parameters  of  the  model,  AV,  Vu  and  yro  influence  the 
magnitude  and  changes  of  both  Vpop_open  (volume  increase  due  to  alveolar  recruitment) 
and  Vsat  (total  volume  of  saturated  (fully-distended)  elements).  The  shape  of  the  clinically- 
measured  p-V  curve  represent  the  change  in  V3at.  The  inflection  pressure,  p0)  is  not  only  the 
pressure  at  which  the  local  compliance,  dV/dp,  is  maximum,  but  also  the  pressure  location 
for  a  maximum  rate  of  increase  in  Vpop_open  as  well  as  in  V3at  (Eqs.(8),(10),  Figs.11,12). 

The  mechanistic  model  of  a  TRS  element  presented  in  this  report  consists  of  a  sim¬ 
ple  piston-spring-cylinder  system  with  the  critical  pop-open  pressure  of  the  element  as 
distribution  parameter.  The  pop-open  volume  (=  volume  that  pops  open  at  the  critical 
pressure)  as  well  as  the  spring  constant  are  assumed  constant  and  common  to  all  elements. 
More  comprehensive  and  detailed  analyses  of  clinical  data  as  well  as  advice  from  clinical 
experts  are  needed  to  advance  the  model  further  and  also  to  make  it  a  practical  tool  for 
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understanding  various  respiratory  conditions.  However,  it  is  believed  that  the  analyses 
presented  here  show  the  developments  and  use  of  a  mechanistic  model  as  a  possible  new 
approach  to  investigate  respiratory  systems. 


Figure/Table  Captions 


Table  1.  Summary  of  Deflation  Data  Analysis. 

Table  2.  Parameters  of  Inflation  Data  Sets  Examined 

Fig.  1.  Data  points  for  inflation  (unfilled)  and  deflation  (filled),  and  the  corresponding 
p-V  equation  of  mechanistic  model  for  inflation  (I).  and  the  error  function 
p-V  equation  for  deflation  (D), 

Pid  —  pressure  at  the  intersect  of  the  inflation  curve  and  a  line  parallel  to  the  x-axis 
passing  through  the  initial  deflation  data  point. 

(See  Table  1  of  Part  I  (II)  for  numerical  values  of  parameters  for  inflation  (deflation).) 

Fig.  2.  Vy  (predicted  from  the' mechanistic  model)  vs  V{f  (of  error-function  p-V  equation 
for  deflation). 

Fig.  3.  VJy  (predicted  from  the  mechanistic  model  with  yro  =  0)  vs  Vjj  (of  error-function 
p-V  equation  for  deflation). 

Fig.  4.  AVd  (predicted  from  the  mechanistic  model)  vs  AVd  (of  error-function  p-V 
equation  for  deflation).  Letters  in  the  figure  =  Data  No. 

Fig.  5.  Deflation  curves. 

triangle  =  data  points,  dotted  =  error  function  p-V  equation  with  the  method  of 
least  squares  applied  to  determine  ( Ad ,  pd,  AVd,  V$), 

solid  =  error  function  p-V  equation  with  (AVd,  V$)  determined  from  Eqs.(4,5) 
and  ( Ad ,  p$)  determined  by  the  method  of  least  squares. 

Fig.  6.  Ranges  of  parameters  of  inflation  data  sets. 

(a)  po  vs  A,  (b)  AV  vs  A,  (c)  1/(1  +  pro)  vs  A.  Letters  in  the  figure  =  Data  No. 

Fig.  7.  Model-based  p-V  equation. 

(a)  Data  Set  B,  E  and  N.  (b)  Data  Set  E  and  T. 
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(c)  Data  Set  N  and  R.  (d)Data  Set  E  and  M. 

Vertical  lines:  dotted  =  po,  broken  =  po  Pro- 

Fig.  8.  Non-dimensional  ( p  —  V )  equation. 

(a)  Data  Set  B,  E  and  N.  (b)  Data  Set  E  and  T. 

(c)  Data  Set  N  and  R.  (d)Data  Set  E  and  M. 

Fig.  9.  Distribution  of  elements. 

(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)). 

Fig. 10.  Number  fraction  of  total  open  elements,  Nopentotai  (p  =  p)/N,  vs  pressure  (solid) 
and  number  fraction  of  saturated  elements,  Nsat_  (p  =  p)/N,  vs  pressure  (broken). 

(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)). 

Fig. 11.  Number  fraction  of  open  elements  at  p  =  0,  Nopen(p  =  0 )/N  vs  A 

Fig. 12.  p  -  V  curve  (solid,  Eq.(2)),  p  -  Vpop-open  curve  (dotted,  Eq.(8a)) 
and  p  —  Vsat  curve  (broken,  Eq.(8b))  for  six  data  sets. 

(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)). 

Letters  in  the  figure  =  Data  No. 

Fig.13.  dVpop_open/dp  [. L/cmH20 ]  vs  p  [cmH20]  (solid) 

and  dV sat/ dp  [L/cmH20]  vs  p  [cmH20]  (broken)  for  six  data  sets. 

(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)). 
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Table  1.  Summary  of  Deflation  Data  Analysis 


Data  No 

Po 

[cmH^O] 

AVd 

[L] 

vi 

[L] 

VS 

[L] 

A  1. 

1.4712 

11.410 

1.9214 

-0.1528 

1.7685 

2. 

1.7567 

1.7064 

B  1. 

1.5528 

12.975 

1.4115 

-0.1227 

1.2887 

2. 

1.3405 

1.2808 

C  1. 

2.3293 

14.808 

1.1753 

0.0800 

1.2599 

2. 

1.2525 

1.2259 

D  1. 

1.0772 

6.846 

2.7012 

-0.6096 

2.0915 

2. 

2.5850 

2.0572 

E  1. 

,  2.2044 

-15.602 

0.7009 

0.0703 

0.7712 

2. 

0.7784 

0.7329 

F  1. 

1.7995 

14.418 

1.0476 

-0.0527 

0.9948 

2. 

1.0041 

0.9252 

G  1. 

0.9839 

8.140 

1.2811 

0.0454 

1.3265 

2. 

1.6496 

1.3455 

H  1. 

0.7147 

"6.472 

2.0772 

-1.0349 

1.0423 

2. 

1.1910 

1.0053 

I  1. 

2.1010 

14.157 

1.4798 

0.5637 

2.0435 

2. 

1.9943 

1.9784 

J  1. 

1.6600 

13.404 

1.8960 

-0.1250 

1.7103 

2. 

1.7759 

1.6858 

K  1.  1 

1.0577 

9.778 

1.5070 

-0.3657 

1.1412 

2. 

1.1620 

1.0860 

L  1. 

0.3107 

3.151 

1.3930 

-0.8076 

0.5853 

2. 

0.8574 

0.6206 

M  1. 

2.1702 

14.989 

1.4322 

0.1442 

1.5764  ■ 

2. 

1.4921 

1.5163 

N  1. 

0.7220 

6.45 

2.4020 

-1.1900 

1.2140 

2. 

1.1654 

1.1115 

0  1. 

1.0558 

10.969 

1.8684 

-0.3195 

1.5488 

2. 

1.6515 

1.4375 

P  1. 

1.5406 

12.231 

0.7461 

0.0684 

0.8145 

2. 

0.8679 

0.8154 

Q  1. 

1.1691 

10.670 

1.1801 

-0.1515 

1.0285 

2. 

1.0938 

0.9857 

R  1. 

0.8643 

8.1411 

1.8260 

-0.4893 

L3366 

2. 

1.7784 

1.2923 

S  1. 

2.7458 

17.050 

2.0015 

-0.5673 

1.4341 

2. 

1.4461 

1.3962 

T  1. 

4.5056 

23.619 

1.1748 

0.1057 

1.2805 

2. 

1.2115 

1.2026 

U  1. 

2.9735 

18.599 

2.3404 

0.0408 

2.3812 

2. 

2.5197 

2.3305 

1.  Obtained  by  method  of  least  squares  with  error  function  p-V  equation 
for  deflation  process. 

2.  Results  from  the  mechanistic  model  of  inflation  process. 
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Table  2.  Parameters  of  Inflation  Data  Sets  Examined 


Data  A 

Po 

AV 

1/(1  +  Vto ) 

Vto 

&D 

Vopen  —  sat.  ( V  —  0) 

B 

2.7304 

21.999 

1.5559 

0.736 

0.359 

12.857 

0.01789 

E 

2.6497 

30.817 

1.6847 

0.676 

0.480 

18.559 

0.02645 

M 

2.9972 

30.327 

4.2463 

0.615 

0.626 

16.147 

0.04932 

N 

2.8046 

15.297 

1.6219 

0.736 

0.358 

8.704 

0.01685 

R 

1.6209 

16.352 

1.7396 

0.711 

0.406 

16.098 

0.07779 

T 

5.4708 

30.037 

1.7694 

0.766 

0.289 

8.761 

0.00012 

Po  and  <7£>  in  [cmH^O] 

7  &V  and  Vopen—sat.  (p 

=  0)  in  [L] 
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V$  (of  error  function  equation  for  deflation) 


Fig  5.  triangle  =  data  points,  dotted  error  function  p-V  equation  with  the  method  of 
least  squares  applied  to  determine  (  Ad  ,pd  ,AVd  ,V„  ) , 

Solid  =  error  function  p-V  equation  with  ( AVd  ,Vd )  determined  from  Eqs.  (4,5)  and  (A d,pd) 
determined  by  method  of  least  squares. 
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Fig.  9.  Distribution  of  elements. 
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(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)) 
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Fig.10.  Number  fraction  of  total  open  elements  N  /  '  Wft, 

v  elements,  l\apen  total  (p  =  p)/jV,  vs  pressure  (sol: 

and  number  fraction  of  saturated  elements  /V  (r,  ^/at 

ents,  iv3at.  [p=p)/N ,  vs  pressure  (broken), 
(a)  through  (f)  for  Data  Seta  B(  =(a)),  E,  M,  N,  R  aud  T(  =(f)). 
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Fig. 11.  Number  fraction  of  open  elements  at  p  =  0,  Nopen(p  =  0)/. 


j)  [cmHiO] 


p  [cmHjO] 

Fig. 12.  p  —  V  curve  (solid,  Eq.(2)),  p  —  Vpop-ovtn  curve  (dotted,  Eq.(8a)) 
and  p  —  Vaat  curve  (broken,  Eq.(8b))  for  six  data  sets. 

(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)). 
Letters  in  the  figure  =  Data  No.  qQ 
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Fig’.13.  dVpop-opcn/dp  [L/cmH20]  vs  p  [ cmH20 ]  (solid) 

and  dV sat/ dp  [L/cmH20]  vs  p  [ cmH20 ]  (broken)  for  six  data  sets, 
(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)).  . 
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Chapter  4.  A  Mechanistic  Model:  Part  III.  Deflation  Process 

Abstract 

The  deflation  process  is  different  structurally  from  the  preceding  inflation  process  because 
of  the  absence  of  the  pop-open  mechanism  and  the  possibility  of  airway  closure  and  collapse. 
An  analysis  is  made  on  the  ideal  deflation  process,  Vd-inf .  (p),  which  results  from  the  elastic 
decrease  of  the  volume  starting  from  the  end-of-inflation  pressure.  Comparisons  between 
the  actual  deflation  curve  and  the  corresponding  Vd-inf.  curve  are  made  to  examine  a 
degree  of  airway  closure  and  collapse  in  terms  of  the  parameters  of  the  mechanistic  model. 
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V 


A 

B(P ) 

C(P) 

D(p) 

F(f  =  PlD ) 

k 

kd 

k 

N 

Nj 

P 

P 

Pcj 

Pcj 

PlD 

PlD 

Po 

pi 

pL 

A  p 
V 

Vd—inf.  ip) 
VL{U) 


Nomenclature 

piston  surface  area  on  which  pressure  is  acting 

Pid  +  1  —  Vt) 

2  +  Pid  +  (p~  Pid) Ik ,  Eq.(7). 

P  +  1  ~  Vt)  +  k  •  (1  +  vto),  Eq.(10). 

P  +  1  -  Vto,  Eq.(9). 

inflation  distribution  functions  (Eq.(2))  at  p  =  pID 
spring  constant  [N / m]  for  inflation 
spring  constant  [N/m]  for  deflation 
=  kd/k. 

total  number  of  TRS  elements 
number  of  elements  at  energy  level  j 
pressure  (interpleural  pressure  difference) 

non-dimensional’ pressure,  p/po  —  1 

critical  pressure  at  which  an  element,  j,  ‘pops  open’. 

Pcj  /  PO 

pressure  at  the  intersect  of  inflation  and  deflation  processes 
after  horizontal  shift  of  the  deflation  curve. 

Pid/po  ~  1 

pressure  at  the  inflection  point  for  inflation  (=  (k/As)  yT0) 

pressure  at  the  inflection  point  for  deflation 

pressure  at  the  inflection  point  after  horizontal  shift 

Poriginai  ~  Pa f ter  shift  f°r  horizontal  shift  of  deflation  p-V  curve 

volume 

the  p-V  deflation  equation  due  to  elastic  displacement  only,  Eq.(13) 
lower  (upper)  bound  of  volume  for  inflation 
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V£(U)  lower  (upper)  bound  of  volume  for  deflation 

V deer  ease  volume  decrease  from  the  end  of  inflation  due  to  elastic  displacement 

Vj  volume  of  an  element  j 

Vv-VL  =  NVoiyro  +  1) 

AVd  V$  -  V£ 

AVd  elementalvolumedecreasefromtheendof  inflation,  Eq.(  3) 

Vo  ‘pop-open’  volume  (=  Asyo) 

yj  piston  displacement  of  an  element  j  in  deflation 

Vj  id  displacement  from  inflation-saturated  position  at  he  end  of  inflation  process 

piston  stroke  limit  for  deflation 
Vo  =  Vo /As  ■ 

yr  piston  stroke  limit 

yro  Vt/Vq 

Greek  symbols: 

a  constant  of  proportionality  in  p-V  equation  for  deflation 

A  ql'PqAV  (non-dimensional  parameter) 

Ad  A  fpr  deflation 

A^  A  for  deflation  after  horizontal  shift 

a  (8/7t)1/2/A,  Standard  deviation  (Eq.(3b)) 

Superscript: 

—  related  to  a  single  TRS  element 

d  deflation  process 

Acronym: 

ARDS  acute  respiratory  distress  syndrome 

TRS  total  respiratory  system 
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Introduction 

The  data  sets  of  ARDS  patients  (Data  Source  A  [1])  previously  analyzed  for  the 
inflation  process  are  the  only  ones  with  complete  inflation-deflation  p-V  measurements. 
However,  the  inflation  and  the  deflation  curves  were  obtained  separately,  requiring  a  hor¬ 
izontal  shift  of  the  deflation  curve  so  that  the  end-of  inflation  pressure  becomes  equal  to 
the  initial  deflation  pressure.  In  the  previous  discussion,  this  has  been  shown  to  be  an 
effective  method  of  analyzing  the  inflation-deflation  process  as  a  whole.  Changes  in  the 
parameters  of  the  p-V  error-function  model  equation  due  to  the  shift  will  be  presented  first 
in  the  analysis  section  to  follow.  Based  on  the  analyses  of  the  inflation  process,  it  may  be 
seen  that  characteristics  of  the  corresponding  deflation  process  depends  on  the  magnitude 
of  Pid  (=  end-of  inflation  pressure  =  initial  deflation  pressure),  making  it  difficult  to  make 
quantitative  comparisons  between  two  different  deflation  p-V  curves.  (The  data  sets  of 
ARDS  patients  range  in  the  magnitude  of  pid  from  18  to  38  cmH20.) 

The  analysis  based  on  the  mechanistic  model  will  examine  an  ideal  deflation  p-V  curve 
that  may  be  obtained  by  decreasing  the  pressure  from  p/£>.  Since  the  pop-open  mechanism 
is  absent  for  the  deflation  process,  the  p-V  curve  thus  predicted  by  the  model  is  a  vol¬ 
ume  decrease  based  on  elastic  contraction  of  the  elemental  volume.  The  actual  deflation 
curve  and  the  corrrespondign  model-based  p-V  curve  will  be  compared  and  discussed  after 
Analysis. 

Analysis 

The  quasi-static  pressure- volume  curves  are  often  obtained  separately  for  the  inflation 
and  the  deflation  processes.  In  order  to  examine  the  two  processes  simultaneously,  it  is 
necessary  to  postulate  a  certain  relationship  between  them.  Particularly  we  are  interested 
in  the  shape  of  a  deflation  curve  directly  preceded  by  an  inflation  curve.  Therefore,  we 
propose  that  the  shape  of  the  deflation  curve  (measured  separately  from  the  corresponding 
inflation  curve)  remains  the  same  as  it  is  shifted  horizontally  along  the  x  (pressure)-axis 
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until  the  starting  data-point  of  the  deflation  curve  lies  on  the  corresponding  inflation  curve. 
Quantitaively  this  means  that  the  upper  and  lower  volume  asymptotes  of  the  deflation 
process,  V$  and  V*,  remain  unchanged  before  and  after  the  translation  of  the  deflation 
curve.  (A  validity  of  this  assumption  has  been  presented  before.) 

Then,  the  deflation  p-V  equation  before  the  horizontal  translation 


v(p  -  M  ~  (M  ±  YlM 

AVd/2 


becomes,  after  the  translation, 


V(p)  ~  (JYu  ±  YlM 

AVd/2 


' p  —  A p 

■  VL 


where  A p  =  p original  -  Pa, f ter  shift ,  with  the  subscript  “a”  indicating  “after  shift”.  The 
identitity  of  the  right  hand  side  of  the  equations  requires 


Pq  P0a ’ 


Therefore,  the  parameters,  A d  and  p0,  before  and  after  the  horizontal  shift  are  related  by 
the  following  equations: 


As  shown  in  the  inflation  analysis  there  are  two  regions  in  the  inflation  p-V  curves; 
an  upper  pressure  region  where  some  elements  are  saturated  (their  piston  stroke  reached 
its  limit)  and  a  lower  pressure  region  where  all  elements  are  active  with  the  boundary 
pressure  between  the  two  regions,  Pboundary,  being  equal  to  pQyT0  (i.e.  pboundary  =  yT o  -  !)• 
Since  the  magnitude  of  uto  is  less  than  1,  it  is  reasonable  to  impose  a  condition  that  the 
end-of-inflation  pressure,  pID,  is  greater  than  Pboundary ;  that  is,  we  analyze  the  case  of 
Pid  —  Pboundary  Pid  =  the  non-dimensional  end-of-inflation  pressure  =  the  non- 
dimensional  initial  deflation  pressure  =pid/Po  —  !• 
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Compared  to  the  inflation  process,  the  volume  after  deflation  to  p  =  0  often  deviates 
substantially  from  zero.  Since  the  volume  change  is  obtained  by  measuring  the  exhaled 
gas  volume  combined  with  the  end-of  inflation- volume,  the  non-zero  value  of  the  volume 
at  p  =  0  indicates  differences  between  the  deflation  and  the  inflation  process. 

In  an  ideal  deflation  process  of  our  mechanistic  TRS  model  for  the  defla¬ 
tion,  each  element  would  reduce  its  (open)  volume  against  the  force  from  the 
spring  as  the  pressure  decreases.  This  deflation  process  is  analyzed  below. 

The  number  distribution  of  elements  for  a  deflation  process  is  fixed  in  terms  of  the  inflation 
distribution  function  at  the  the  end  of  the  inflation  process  as, 


(2) 


where  Pcj-Pcj/Po,  cr  =  (8/rc)  a  /A . 

For  each  element  its  volume  change  is  to  be  evaluated  as  pressure  decreases,  and  a  sum 
(an  integral)  of  the  elemental  volume  change  over  all  elements  using  the  elemental  number 
distribution  at  the  end  of  the  inflation  process,  Eq.(2),  yields  the  deflation  p-V  equation 
as  the  volume  decreases  from  the  end  of  the  inflation  pressure,  pro  ■ 

During  the  inflation  process,  a  force  balance  for  an  element  j  at  p  =  p  >  pCj  is, 

Asijp  -  pCj)  =  kyj  <  kyr-  (See  Fig.l.)  Therefore,  for  the  element  j  with  its  critical 
pop-open  pressure,  pCj ,  the  piston  stroke  limit  is  reached  when  p  =  ( k/Aa)yT  +  Pcj 
(or  p  =  PoPto  +  Pcj >  recognizing  the  relation,  ( k/A$)y0  =  p0,  obtained  in  the  inflation 
analysis). 

Then,  an  element  at  the  end  of  the  inflation  process  belongs  to  either  one  of  the 
following  two  cases: 

(Case  1)  At  the  end  of  inflation  process  of  p  =  pm ,  elements  are  unsaturated  (piston 
remaining  active)  if  their  critical  pop-open  pressures  are  in  the  range  of 

PID  ~  Po  Pro  <  Pcj  <  P/D- 


107 


(Case  2)  An  element  is  inflation-saturated  at  p  —  pid  if  its  critical  pop-open  pressure, 
pCj ,  is  in  the  range  of 

0  <  pCj  <  Pid  ~  PoVto- 

Therefore,  the  volume  change  during  the  ensuing  deflation  process  must  be  evaluated  for 
the  two  cases  separately. 

Case  1:  p/d  -  Po  pro  <  Pcj  <  Pid  (he.  A  <pCj  <  Pid  +  1  with  A  =  pID  +  1  -  pro) 

(Elements  that  are  not  inflation-saturated  at  p/d) 

Referring  to  Fig. 2  the  force  balance  at  p  =  p  (<  pID)  and  at  p  —  pID  are  respectively, 

Pid~P  =  -^(Vj  -Vjw)>  Pid  -  Pcj  =  j-  (wr  -  vfm) 

where  kd  =  spring  constant  for  deflation,  yjID  =  the  displacement  from  the  inflation- 
saturated  position  at  the  end  of  the  inflation  process  of  p  =  p/D-  They  yield  the  following 
equations  for  AVf  (=  an  elemental  volume  change  (decrease)1  from  the  end  of  the  inflation 
process  =  A3(fff  -  yjID))  in  terms  of  pressure  change  from  the  end  of  the  inflation  as  well 
as  Pj  in  terms  of  the  inflation  pop-open  pressure,  pCj) 

AVf=A(pw-p),  fl  =  A  (3,4) 

The  range  of  yf  is  limited  by  the  imposed  condition  of  yj  <  (=  piston  stroke  limit 

for  deflation).  Therefore,  Eq.(4)  yields  the  condition  for  element  j  to  reach  the  state  of 
deflation-saturation  at  p  =  p  as 

Pci  =  T,(pi~:ST  +  ^p~MTd~  l)pw )  •  (5) 

Also,  after  the  stroke  limit  is  reached  the  elemental  volume  change  stops  at 

AVf  =  As(pf  -  pr  +  -~(pid  ~  Pcj))-  (6) 

In  summary,  for  the  elements  not  inflation-saturated  (pid  —  PoVto  <  Pcj  <  Pid)  the 
stroke  limit  of  yf  is  reached  (deflation-saturated)  for  elements  satisfying  Eq.(5);  while, 
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AV^  for  deflation-unsaturated  elements,  and  deflation-saturated  elements  are  respectively 
given  by  Eqs.(3)  &  (6). 

In  non-dimensional  representation  and  referring  to  Fig. 3,  for  elements  with 
A  <  pcj  <  PID  +  1, 

(l-A)If  pcj  >  B(=2  +  pID  +  (p-pID)/k,  k  =  kd/k), 
the  element,  j,  is  deflation-saturated  at  p  —  p 

with  AVf  =  V0  (2  +  p!D  -  pcj). 

(1-B)  If  %j  <  B, 

the  element,  j,  is  deflation-unsaturated  at  p  =  p 

-  with  AVf  =  V0(p-  pID)/k. 

Case  2:  '  0  <  pcj  <  pID  -  p0  yT0  (i.e.  0  <  pcj  <  A) 

(Elements  that  are  inflation-saturated  at  pm) 

As  shown  in  the  inflation  analysis  an  open  element,  j,  becomes  inflation-saturatred 
when  pressure  during  the  inflation  reaches  (pCj  +  PoVto),  which  is  also  the  pressure  at 
which  the  element  becomes  free  dining  the  deflation  (deflation-unsaturated).  (See  Fig.4 
(a),  (b).)  Therefore,  at  p  =  p  in  the  deflation  process, 

the  element,  j,  is  inflation-saturated  if  0  <  pcj  <  p  —  po  yro- 
the  element,  j,  is  deflation-unsaturated  if  p  —  po  Vto  <  Pcj  • 

From  the  force  balance  at  p  =  p  for  a  deflation-unsaturated  element,  the  elemental  volume 
change  during  the  process  from  pin  to  p  is  shown  to  be 

AVf  ~  Agfij 

=  (A2Jkd)  (Pcj  +  Po  Vto  -  p)- 

The  element  becomes  deflation-saturated  at  y*f  =  (Fig.4  (c)).  A  force  balance  applied 
to -an  element  at  the  onset  of  the  deflation  saturation  is 

(Pc  j  +  Po  Vto)  -P-  ( kd/As )  •  y%. 


(7) 


(8) 
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Hence,  the  elements  with  pcj  >  ( kd/As )  -yf+p-po  yro  are  deflation-saturated  with  its 
volume  change  being  stopped  at  AVd  =  Vo  (1  +  yro)- 

In  summary,  for  an  element,  j,  with  their  critical  pop-open  pressure  in  the  range  of 
0  <  Pcj  <  A  (element  that  is  inflation-saturated  at  Pw) 

(2- A)  If  0  <  pCj  <  D(p)  (=p  +  1  -  yro), 

the  element  remains  inflation-saturated  at  p  =  p 

with  A  Vf  =  0,  (9) 

(2-B)  If  D  <  pCj  <  C(p)  (=  p  +  1  —  yro  +  k  •  (1  +  yro)) 
the  element  is  deflation-unsaturated  at  p  =  p 

VS 

with  A  Vf  =  5  (pCj  +  yro  ~(p  +  1)),  (10) 

k 


(2-C)  If  C  <  pcj  <  A, 

the  element  is  deflation-saturated  at  p  =  p 


with  AV/  =  V0(l+yT0), 


(11) 


A  sum  of  the  elemental  volume  change,  A Vd,  over  all  elements  yields  the  volume 
decrease  from  the  end  of  the  inflation  process  due  to  elastic  displacement  only,  Vfecrease(p)] 
that  is, 


vd 

v  r 


decrease 


(p)  =  J  A  VfdNj 


all  Nj 


=  J  N  AVf  F(pID)  dpCj 


all  Pc  j 
rPiD+l 

Jo 

AV 
1  +  yro 


-  /AVd\ 

NV0[-^-)F(pID)dp 


"A  -AVd 


rA  /  AV-\ 

/  \~lTJn  /^d)  dPcj 

JO  '  Vq  '  Case 2 


+  r,D+1  (S) 

Ja  V  Vn  ' 


Yq  '  Case  1 


F(pID)dp 


CJ 


(12) 
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where  (AVf/Vo)  may  be  found  from  Eqs.(7)-(11)  as  function  of  pcj  and  p. 

Then,  Vd-inf(p ),  defined  as  the  p-V  deflation  equation  due  to  elastic  displacement 
only,  may  be  expressed  as 

Vd-inf{p )  =  VinflationiP  —  Pw)  ~  V decrease (P )  (13) 

where  Vinfiation(p  =  Pjp)  is  the  volume  at  the  end  of  the  preceding  inflation  process. 
Because  of  the  complexity  associated  with  the  integrals  in  Eq.(12),  Vd-inf(p)  is  evaluated 
computationally.  Also,  the  ratio  of  fspring  constant,  k,  between  the  inflation  and  the 
deflation  process  needs  to  be  specified.  Since  the  idealized  deflation  process,  Vd-inf  vs  p 
curve,  starts  at  the  end-of-inflation  pressure,  and  ends  at  p  =  0,  with  Vinjiation  =  Vd-inf 
at  both  ends  of  the  process,  the_magnitude  of  k  is  determined  by  satisfying  the  condition 

of  Vinflation{p  =  0)  =  Vd—inf(p  =  0). 

The  computational  procedure  is  outlined  below. 

Input  Data:  A,  AV,  p0,  VL,  pID,  yT0. 

1.  Find  the  value  of  &. 

1-1.  Assign  an  initial  guess  value  for  k. 

1-2.  Evaluate  Vd-inf(p  =  0)  from  Eqs.(12),(13). 

1-3.  If  |  Vd-inf{jp  =  0)  -  Vinfiation{p  =  0)  |  <  e  (=  convergence  criterion  <<  1), 
k  is  found. 

If  not, 

let  k(new )  =  k(old)  +  Ak  (A k  =  increment  for  computational  iteration) 

and 

repeat  Step  1-2  and  1-3  until  convergence  is  achieved. 

2.  Evaluate  p  -  Vd-inf  curve  from  Eqs.(12),(13)  over  the  range  of  0  <  p  <  pID. 
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Results  and  Discussion 

As  noted  before,  the  volume  change  in  the  deflation  process  is  obtained  by  subtracting 
the  measured  exhaled  gas  volume  from  the  end-of-inflation  volume;  hence,  the  non-zero 
value  of  the  volume  at  the  end  of  the  deflation  process  at  p  =  0  may  be  indicative  of  the 
gas  volume  trapped  in  the  TRS  caused  by  such  effect  as  airway  closure.  We  evaluate  the 
volume  change,  Vd~inf,  of  an  idealized  deflation  process  in  which  the  elemental  volume 
change  from  the  conditions  at  the  end  of  the  inflation  process  is  solely  due  to  the  elastic 
deflation  (=  piston  displacement  as  the  spring  expands  with  a  decrease  in  pressure). 

Table  1  lists  inflation  parameters  of  the  mechanistic  model  (A,  p0,  AV,  yro),  AVd 
(=  AV  of  deflation),  A^,  and  p^a  (A  and  po  of  the  deflation  process  after  the  horizontal 
shift),  a  ratio  of  the  volume  at_the  end-of-inflation  pressure,  V (pjd),  to  AVd,  and  the 
spring  constant  ratio,  k  (  kd/k  —  k deflation /^inflation )  foe  twenty  one  sets  of  Data 
Source  A.  The  p-V  deflation  curves  of  the  twenty  one  sets  are  also  shown  in  Appendix, 
in  which  the  solid  curve  is  the  deflation  curve  after  the  measured  p-V  deflation  curve  is 
horizontally  translated  so  that  Pend— of —inflation  ~  Pinitialdef lation  ]  while,  the  dotted  curve 
is  the  computed  results  of  the  ideal  p-V  deflation  curve  predicted  from  the  mechanistic 
model. 

Both  A  and  p0  decrease  from  the  inflation  to  the  deflation  process  as  a  result  of  the 
structural  differences  between  the  two  processes.  Negative  values  for  Ad,  are  shown  in 
some  data  sets  of  Table  1  due  to  a  horizontal  shift  of  the  original  deflation  curve  to  the 
left  on  p-V  diagram.  (See  Eq.(l).) 

Data  G  in  Table  1  is, 

Data  A  po  AV  yT0  AVd  A*  pda  V(pID)/AVd  k 
G.  1.9277  20.731  2.8704  0.447  1.281  -0.4288  -3.547  1.008  0.30 

Fig. 5  shows  Data  G  of  Table  1  (the  p-V  deflation  curve  after  the  horizontal  shift  (solid)  and 
the  corresponding  ideal  p-V  deflation  curve  (dotted)).  A  large  positive  value  of  Vdef{p  =  0) 
suggests  a  substantial  airway  closure  during  the  deflation.  In  terms  of  numerical  values  it 
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Other  data  sets  with  negative  pdQa  are  Data  H,  L,  N,  shown  in  Fig.7,  all  characterized 
by  a  small  pressure  range  compared  to  the  data  with  positive  p$a.  The  deflation  curve 
(after  the  horizontal  shift)  of  Data  L,  in  particular,  indicates  small  volume  change  with 
the  decrease  in  pressure;  which  is  due  to  a  large  negative  value  of  p$a.  Compared  to  Data  G 
their  values  ofV(pID)/AVd  are  smaller  than  1.0.  Also  noted  is  an  increase  in  magnitude 
from  AV"  to  AVd  for  all  three  sets  with  a  negative  value  of  pda . 


Data 

A 

Po 

Ay 

Vto 

AVd 

AJ 

Poa 

V(pID)/AVd 

k 

H. 

2.0421 

15.405 

1.7129 

0.431  , 

2.077 

-0.6247 

-5.657 

0.466 

0.30 

L. 

1.5318 

13.213 

1.1256 

0.695 

1.393 

-1.8841 

-18.696 

0.411 

0.41 

N. 

2.8046 

15.297 
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Fig.  7  Deflation  analysis  of  Data  H,  L,  N 
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Data  D  and  K  are  similar;  i.e.  a  small  change  between  AV  and  AVd  with  the  deflation 
parameters,  AQ  and  Pga,  being  roughly  the  same.  However,  the  difference  between  the  Vd 
and  Vd_inf  is  large  in  Data  D  and  very  small  in  Data  K.  In  other  words,  the  difference 
between  two  curves  in  D  and  K  is  due  to  a  magnitude  of  k,  which  affects  the  shape  of 

the  p -Vd-inf  curve.  A  high  value  of  k  in  Data  K  results  in  negligible  effects  of  the  airway 
closure. 
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AV 
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AVd 

Ad 

.a 

POa 

V(Pid)/ AVd 

k 

D. 

1.9257 

13.999 

2.8392 

0.379 

2.701 

0.7625 

4.846 

0.752 

0.27 

K. 

2.4708 

17.951 

1.3304 

0.440 
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0.6000 

5.549 

0.645 
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Fig.  8  Deflation  analysis  of  Data  D,  K 
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Finally  discussed  is  Data  T,  a  data  set  with  very  high  values  for  both  A  and  po  in  the 
inflation  curve  as  well  as  for  both  Ad  and  pda  in  the  corresponding  deflation  curve.  This 
is  the  only  data  in  which  the  solid  p  —  V^after  horizontal  shift)  curve  is  located  below  the 
dotted  ideal  p  —  V  curve  over  a  majority  of  the  pressure  range.  This  may  imply  that  the 
airway  as  well  as  alveoli  are  collapsing  at  a  rate  faster  than  the  elemental  elastic  deflation 
predicted  by  the  ideal  curve  with  a  very  small  amount  of  trapped  air. 

Data  A  p0  AV  yT0  AVd  A*  pda  V(pw)/AVd  k 
T.  5.4708  30.037  1.7694  0.289  1.175  3.9600  20.759  1.037  0.22 
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Summary  and  Conclusions 

•  An  analysis  is  made  on  the  ideal  deflation  process,  Vd-inf.  (p) ,  which  results  from  the 
elastic  decrease  of  the  volume  starting  from  the  end-of-inflation  pressure  in  the  mechanistic 
model.  A  comparison  between  the  actual  deflation  curve  and  the  corresponding  Vd-inf. 
curve  is  made  to  examine  a  degree  of  airway  closure  and  collapse  in  terms  of  the  parameters 
of  the  mechanistic  model.  A  very  large  value  of  Ad  (>  4)  is  accompanied  by  large  airway 
collapse  without  trapped  air;  while,  the  data  with  k  (ratio  of  spring  constant  between  the 
deflation  and  the  preceding  inflation)  of  the  order  of  unity  shows  negligible  effects  of  airway 
closure  on  the  deflation  curve.  The  differences  between  the  measured  deflation  curve  and 
the  corresponding  ideal  deflation  curve,  Vd-inf. ,  evaluated  from  the  mechanistic  model  as 
a  degree  of  airway  closure  depend  on  the  relative  magnitude  of  AV  and  AVd,  of  pf)a  and 
P/d,  as  well  as  the  magnitude  (positive  or  negative)  of  pda. 

The  deflation  process  yields  a  new  set  of  information  that  may  not  be  obtained  from  the 
analyses  of  the  inflation  process.  Results  obtained  in  this  chapter,  if  combined  with  those 
from  previous  chapters,  may  distinguish  p-V  data  . sets  beyond  what  may  be  achieved  by 
examinations  of  their  original  raw  data.  In  this  chapter  our  data  analyses  are  limited  to  the 
data  sets  from  ARDS  patients  [1]  because  of  data  availability.  The  p-V  curve  analyses  had 
been  made  in  the  past  to  examine  normal  and  diseased  lungs  [2,  5],  patients  with  chronic 
airflow  obstruction  (COPD)  [3],  of  asthma  [4],  as  well  as  an  indicator  for  the  protective 
ventilator  strategy  [6,7].  The  next  step  in  our  study  should  be  aimed  at  establishing  a 
relationship  between  various  respiratory  diseases,  their  p-V  curve  characteristics  and  the 
predictions  based  on  the  mechanistic  model. 
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Table  1.  Summary  of  Data  Analysis. 


Data 

A 

Po 

Ay 

Vto 

Ayd 

Ad 

POa 

V(pw)/AVd 

k 

A. 

2.9578 

22.411 

2.3772 

0.347 

1.921 

1.0899 

8.453 

0.881 

0.26 

B. 

2.7304 

21.999 

1.5559 

0.359 

1.412 

1.1220 

9.375 

0.889 

0.52 

C. 

3.3664 

25.073 

1.6139 

0.329 

1.175 

1.7935 

11.402 

1.055 

0.24 

D. 

1.9257 

13.999 

2.8392 

0.379 

2.701 

0.7625 

4.846 

0.752 

0.27 

E. 

2.6497 

30.817 

1.6847 

0.480 

0.701 

1.8794 

13.302 

1.024 

0.32 

F. 

2.6423 

23.863 

1.4874 

0.365 

1.048 

1.5973 

12.798 

0.863 

0.26 

G. 

1.9277 

20.731 

2.8704 

0.447 

1.281 

-0.4288 

-3.547 

1.008 

0.30 

H. 

2.0421 

15.405 

1.7129 

0.431 

2.077 

-0.6247 

-5.657 

0.466 

0.30 

I. 

3.5449 

25.232 

2.7508 

0.368 

1.480 

2.0066 

13.521 

1.338 

0.27 

J. 

2.7364 

26.830 

3.9129 

0.474  . 

1.896 

1.2389 

10.004 

0.816 

0.34 

K. 

2.4708 

17.951 

1.3304 

0.440 

1.173 

0.6000 

5.549 

0.645 

1.11 

L. 

1.5318 

13.213 

1.1256 

0.695 

1.393 

-1.8841 

-18.696 

0.411 

0.41 

M. 

2.9972 

30.327 

4.2463 

0.626 

1.432 

1.7691 

12.219 

0.987 

0.38 

N. 

2.8046 

15.297 

1.6219 

0.358 

2.402 

-0.4079 

-3.644 

0.458 

0.26 

0. 

2.1412 

18.837 

1.7753 

0.389 

1.868 

1.1520 

11.969 

0.755 

0.28 

P. 

2.4553 

27.032 

1.3278 

0.467 

0.746 

1.3643 

10.831 

1.057 

0.31 

Q. 

2.2041 

19.583 

1.2941 

0.389 

1.180 

0.9500 

8.670 

0.823 

0.28 

R. 

1.6209 

16.352 

1.7396 

0.406 

1.826 

0.9705 

9.141 

0.687 

0.28 

S. 

3.0894 

26.925 

3.1878 

0.407 

2.002 

1.3818 

8.580 

0.672 

0.28 

T. 

5.4708 

30.037 

1.7694 

0.289 

1.175 

3.9600 

20.759 

1.037 

0.22 

U. 

3.2819 

24.349 

2.8956 

0.327 

2.340 

2.8456 

17.799 

0.995 

0.24 

1.  V  and  Ay  in  [L],  P0  in  [cmH20]. 

2.  Inflation  parameters  from  the  mechanistic  model. 
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Deflation  at  p  <  plD  (deflation  saturated)  Deflation  at  p  <  plD  (deflation  unsaturated) 

(l-A)  (1-B) 


Fig.3.  Summary  of  Case  1.  E.O.I  =  End  Of  Inflation 
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P=P„ 


E.O.I  (Element  at  End  of  Infaltion) 


Deflation-unsaturated  element 


Deflation  -saturated  element 


Fig  .4.  Case  2  (Elements  that  are  inflation-saturated  at  p  =  plD) 
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Appendix 


e  p-V  deflation  curve  after  the  horizontal  shift  (solid)  and 
the  corresponding  ideal  p-V  deflation  curve  (dotted)) 
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Key  Research  Accomplishments 

1.  The  representation  of  p-V  curves  by  a  single  (non-linear)  model  equation  (either  the 
error  function  p-V  model  equation  or  the  sigmoidal  model  equation)  is  confirmed  to 
be  an  effective  method  for  clinical  data  analyses. 

2.  A  mechanistic  model  is  constructed  for  the  inflation  process;  which  makes  it  possible 
to  predict  the  internal  conditions  of  the  respiratory  system  from  the  p-V  curve.  The 
model  of  the  inflation  process  also  predicts  the  upper  volume  limit  of  the  deflation, 
Vy,  the  magnitude  of  which  agrees  well  with  the  values  determined  from  the  measured 
deflation  curves;  thus,  justifying  the  method  of  our  model  development. 

3.  A  mechanistic  model  is  derived  for  an  ideal  deflation  process  in  which  the  elastic  (non¬ 
pop-open)  contraction  is  the  only  mode  of  volume  decrease.  The  ideal  p-V  curve  is 
compared  to  the  corresponding  measured  p-V  deflation  curve  with  their  differences 
interpreted  in  terms  of  the  inflation  and  deflation  parameters  in  relation  to  airway 
closure  and  collapse  during  the  deflation  process. 
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Conclusions 


It  is  confirmed  that  both  the  sigmoidal  (tangent  hyperbolic)  p-V  model  equation  and 
the  error  function  model  equation  represent  quasi-static  p-V  curves  extremely  well.  The 
error  function  p-V  model  equation  is  of  the  form, 


V  -VL 
AV 


where  AV  =  V\j  —  Vl,  A  =  ap0  AV(a  =  proportionality  constant),  p  =  p/p0  —  1.  The 
curve  varies  smoothly  between  the  low  pressure  asymptote,  Vl,  and  the  high  pressure 
asymptote,  Vy.  The  midpoint  of  the  curve  is  the  inflection  point  with  its  pressure  denoted 
by  pq.  All  clinical  p-V  curves  made  available  to  us,  ranging  from  healthy  humans  to  ARDS 
patients,  are  represented  well  by  the  model  equation  with  differences  between  data  sets  ex¬ 
pressed  quantitatively  in  terms  of  the  magnitudes  of  the  non-dimensional  parameters  in  the 
model  equation  as  well  as  the  pressure  range  of  the  data  relative  to  the  midpoint  pressure, 
Pq.  Mathematical  processing  of  data  using  the  continuous  model  equation  makes  exami¬ 
nations  of  p-V  curve  characteristics  easier  and  more  accurate,  compared  to  the  piecewise 
representations  by  previous  investigators.  We  found  from  the  literature  survey  and  data 
examinations  that  there  is  no  standardized  procedure  for  p-V  curve  measurements.  For  ex¬ 
ample,  some  published  data  cover  only  a  very  narrow  pressure  range;  while,  others  do  not 
measure  the  esophagus  pressure  as- the  lung  pressure  (not  the  interpleural  pressure)  is  used 
for  the  pressure  scale.  Differences  in  the  p-V  curves  may  be  observed  qualitatively  upon 
visual  comparisons  between  a  healthy  human  and  an  ARDS  patient;  however,  in  terms 
of  the  model  parameters,  quantitative  differences  may  be  detected  among  healthy  adults 
as  well  as  among  ARDS  patients  that  may  not  be  accomplished  by  the  visual  inspection 
alone. 

A  mechanistic  model  of  TRS  (total  respiratory  system)  elements,  each  consisting  of  a 
piston-spring  subsystem  in  a  chamber,  is  developed  to  relate  the  p-V  curve  characteristics 
to  the  internal  change  of  the  corresponding  TRS.  The  mechanistic  model  accommodates 
both  the  alveolar  recruitment  (in  terms  of  the  pop-open  pressure)  and  the  elastic  dis¬ 
tension  of  wall  tissues  (in  terms  of  piston  displacement).  A  (critical)  pressure  at  which 
an  element  (pop-)  opens  is  different  from  element  to  element;  which  yields  the  following 
normal  distribution  for  a  large  number  of  elements  over  the  critical  pop-open  pressure: 


dNj 

N  ■  dpCj 


=  m, 


m  = 


\/2t r  t 


ex 


P(~  of 


1  rPcj  P  -\2 


n- 


where  dNj  =  number  of  elements,  the  critical  pop-open  pressure  of  which  lies  between  pCj 
and  pCj  +  dpCj ,  N  —  total  number  of  TRS  elements  and  pcj  =  pCj/po ■  The  distribution 
has  its  mean  at  p  (=  p/po  —  1)  and  its  standard  devistion,  cr,  being  related  to  the  non- 
dimensional  parameter,  A,  of  the  p-V  model  equation  as  cr  =  (8/7t)s /A.  Each  element  at 
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a  specified  pressure  is  in  a  state  between  closed  conditions  and  fully-distended  (saturated) 
conditions.  As  the  pressure  (p)  increases  along  an  inflation  process  more  elements  are 
recruited  and  open  elements  are  distended,  thus  increasing  the  number  of  open  as  well  as 
saturated  elements.  The  magnitude  of  p0  (pressure  at  the  maximum  compliance)  and  its 
location  relative  to  the  range  of  the  measured  p-V  curve  are  shown  to  be  important  factors 
affecting  changes  of  the  distribution  with  pressure.  Other  parameters  of  the  p-V  model 
equation  influence  the  magnitude  and  changes  of  both  Vpop_open  (volume  increase  due  to 
alveolar  recruitment)  and  Vsaturated  (total  volume  of  saturated  (fully-distended)  elements), 
two  quantities  important  for  the  optimization  of  ventilator  strategy.  Parameters,  V$  and 
AVd  (superscript  indicating  the  deflation  process)  that  are  predicted  from  the  mechanistic 
model  for  the  preceding  inflation  process  agree  well  with  those  of  the  measured  p-V  curves 
of  the  deflation;  thus  providing  justification  for  the  validity  of  the  mechanistic  model. 

The  (pop-open)  alveolar  opening  mechanism  is  absent  for  the  deflation  process;  in¬ 
stead,  the  closure  and  collapse  of  the  airways  affect  the  p-V  curves.  The  modeling  approach 
to  the  deflation  process  is  to  compute  an  ideal  process  of  deflation  by  the  elastic  contrac¬ 
tion  of  TRS  elements  as  the  pressure  is  decreased.  The  difference  in  the  volume  decrease 
between  the  measured  and  the  ideal  deflation  curve  indicates  the  effects  of  the  airway  clo¬ 
sure  on  the  deflation  process.  The  deflation  model  is  computationally  performed  due  to 
complexities  associated  with  the  integrals  involved.  Comparisons  between  the  measured 
deflation  curves  of  ARDS  patients  and  the  corresponding  ideal  deflation  curve,  evaluated 
from  the  mechanistic  model  indicate; 

A  very  large  value  of  Ad  (>  4)  is  accompanied  by  large  airway  collapse  without 
trapped  air. 

A  data  set  with  k  (ratio  of  spring  constant  of  TRS  elements  between  the  deflation 
and  the  preceding  inflation)  of  the  order  of  unity  shows  negligible  effects  of  airway 
closure. 

A  degree  of  airway  closure  depends  on  the  relative  magnitude  of  AV  and  AVd,  of  pd 
and  piDt  as  well  as  the  magnitude  (positive  or  negative)  of  p$. 

It  is  believed  that  the  analyses  presented  in  this  report  show  the  developments  and  use 
of  the  mechanistic  model  as  a  new,  effective  approach  to  investigate  respiratory  systems 
of  patients  with  acute  lung  injury.  Comprehensive  data  analyses  based  on  the  mechanistic 
model  are  made  for  twenty  one  data  sets  of  ARDS  patients,  the  only  available  data  sets 
with  complete  inflation  and  deflation  p-V  curves.  The  p-V  curves  had  been  utilized  in  the 
past  to  examine  normal  and  diseased  lungs,  patients  with  chronic  airflow  obstruction,  of 
asthma,  as  well  as  an  indicator  for  the  protective  ventilator  strategy.  The  next  step  in  our 
study  should  be  aimed  at  establishing  quantitative  and  comprehensive  relations  between 
clinical  diagnoses  of  various  respiratory  diseases  (including  the  underlying  disease  of  ARDS 
patients),  their  p-V  curve  characteristics  and  the  predictions  based  on  the  mechanistic 
model. 
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